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Abstract: 


Acknowledgment : 


Results  of  a study  of  fully  cavitating  hydrofoil  sections  are 
reported.  All  calculations  are  based  upon  the  linearized  theory 
of  cavity  flow  in  two  dimensions.  This  report  is  the  first 
detailed  exploration  of  the  consequences  of  the  general  inverse 
theory  which  permits  the  designer  to  specify  the  design  values 
of  the  lift  coefficient,  cavitation  number  and  the  thickness  of 
the  upper  surface  of  the  cavity  at  the  profile  trailing  edge  as 
well  as  the  shape  of  the  pressure  distribution  on  the  wetted 
surface.  The  ordinates  of  the  upper  cavity  contour  and  the 
wetted  surface  contour  are  calculated.  The  design  angle  of 
attack,  the  cavity  length,  the  drag  coefficient  and  the  moment 
coefficient  are  also  calculated.  It  is  found  for  almost  any 
cavitation  number  and  any  design  lift  coefficient  that  if  the 
center  of  pressure  is  placed  as  closely  as  possible  to  the 
profile  leading  edge  the  resulting  profile  will  have  the  most 
favorable  lift-to-drag  ratio.  The  study  also  includes  off-design 
calculations,  in  accordance  with  the  direct  theory  of  cavity 
flows,  to  determine  cavity  interference  with  the  upper  nonwetted 
surface  of  the  profile  and  the  hydrodynamic  forces  of  particular 
designs. 

This  report  covers  one  of  the  tasks  performed  under  the  NSRDC/ARL 
Cooperative  Hydrofoil  Program  in  FY  1975  for  the  Ships  Perfor- 
mance Department,  NSRDC  under  the  direction  of  Mr.  R.  Wermter 
and  Dr.  Y.  Shen,  The  authors  are  grateful  for  helpful  discussions 
of  this  work  by  Dr.  Shen  and  Professor  J.  W.  Holl.  This  work 
was  sponsored  by  the  Naval  Material  Command  under  Program 
Element  62544N,  Task  Area  ZF43  421  001,  administered  by  the  Naval 
Ship  Research  and  Development  Center  under  contract  N00017-73- 
C-1418. 
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Nomenclature 

cavity  parameter 
lift  coefficient 
drag  coefficient 
moment  coefficient 
cavitation  number 
cavity  length,  & *■  a + 1 
lift-to-drag  ratio 

pressure  amplification  factor  used  in  the  second  design 
method 

nondimens ionalized  perturbation  pressure 
cavity  pressure 
free  stream  pressure 

peak  pressure  location,  measured  from  the  nose 
upper  (non-wetted)  surface  of  the  hydrofoil 
cavity  thickness  at  the  trailing  edge 
free  stream  velocity 

function  used  to  compute  the  angle  of  attack  corresponding 
to  the  cavity-foil  interference  in  the  off-design  calcu- 
lations 

absissa,  taken  as  positive  downstream,  with  the  origin  at 
the  nose 

distance  of  the  center  of  pressure  from  the  nose 
cavity  contour  (upper  surface) 

cavity  thickness  due  to  the  complementary  function  (not 
total  thickness) 
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INTRODUCTION 


General : 


In  1964,  a method  for  calculating  the  shape  of  a fully-cavitating 

hydrofoil  section  with  prescribed  pressure  distribution  was  published 
, * 

[1,  2,  3] . This  work,  which  is  bated  upon  the  linearized  theory  of 
cavity  flews  in  two  dimensions,  resulted  in  three  profile  design  pro- 
cedures. In  all  three,  the  shape  of  tne  pressure  distribution  on  the 
wetted  surface  of  the  profile  is  given,  and  the  cavitation  number , the 
cavity  thickness  at  the  profile  trailing  edge,  and  the  total  section 
lift  coefficient  are  prescribed.  Then  the  shape  of  the  wetted  surface 
is  calculated  along  with  the  cavity  length,  the  drag  coefficient,  and 
the  moment  coefficient  about  the  profile  nose.  The  lift-to-drag  ratio 
then  follows.  The  three  design  procedures  differ  in  the  way  that  the 
cavity  thickness  requirement  is  satisfied.  Because  of  this  difference, 
the  design  attack  angle  is  an  outcome  of  the  calculations  in  the  first 
and  second  procedures.  In  the  third  method,  the  attack  angle  is 
prescribed  from  the  outset. 

In  the  first  design  procedure,  the  cavity  thickness  is  obtained  by 
use  of  the  complementary  function  for  the  problem.  This  complementary 
function  is  commonly  called  the  "point-drag-solufcion”  [4]  of  cavity  flow 
theory.  In  our  procedure,  the  profile  is  taken  to  be  at  the  attack 
angle  for  "ahockless  entry”  so  that  the  prescribed  pressure  distribution 
defines  completely  the  pressure  distribution  on  the  wetted  surface. 

Thus,  profile  design  at  the  ideal  attack  angle  provides  the  designer  with 


Numbers  in  brackets  refer  to  citations  in  the  list  of  references. 
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the  greatest  control  over  every  detail  of  the  pressure  distribution; 
and  for  this  reason  we  might  be  justified  in  regarding  the  first  design 
procedure  as  the  proper  embodiment  of  the  linearized  inverse  hydrofoil 
problem.  The  analogous  method  of  thin  airfoil  theory  is  given  in  [5, 

6,  7],  for  example.  The  fact  that  the  designer  can  now  control  the 
cavity  thickness  allows  him  to  make  sure  that  the  upper  surface  of  the 
cavity  does  not  interfere  with  the  hydrofoil,  and  thereby  to  make  some 
provision  for  the  structural  adequacy  of  his  design  from  the  start. 

The  second  design  procedure  accounts  for  the  fact  that  interference 
between  the  hydrofoil  and  the  cavity  can  also  be  eliminated  by  increasing 
the  attack  angle  of  the  profile  abovo  tne  ideal  angle,  instead  of  by 
means  of  the  complementary  function.  In  this  case,  one  prescribes  the 
value  of  CL  as  before.  However,  now  one  cannot  control  all  details  of 
the  pressure  distribution.  Only  the  shape  of  that  part  of  the  pressure 
distribution  which  is  associated  with  the  wetted-surface  camber  is 
subject  to  the  designer's  pleasure.  The  pressure  distribution  now  has 
two  parts:  the  first  part  which  is  due  to  the  nose  singularity  (or 

incremental  attack  angle),  and  the  second  part  which  i£  due  to  camber 
of  the  wetted  surface.  The  shape  of  the  first  part  is  fixed,  but  we 
can  vary  its  magnitude  by  changing  the  angle  of  attack.  Once  the  shape 
of  the  second  part  is  specified;  one  can  adjust  its  amplitude  so  that 
the  value  of  the  lift  coefficient  originally  specified  is  met  by  the 
sum  of  contributions  frem  the  first  and  the  second  parts  of  the  pressure 
distribution. 


The  third  design  procedure  is  more  general  than  the  first  and 
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second  procedures,  because  the  complementary  function  and  the  incremental 
attack  angle  are  both  used  to  meet  prescribed  cavity-thickness  require- 
ments. Now,  one  specifies  the  value  of  attack  angle  along  with  the  lift 
coefficient,  cavitation  number,  cavity  thickness  at  the  trailing  edge, 
and  the  shape  of  that  part  of  the  pressure  distribution  due  to  camber. 
Then,  the  remaining  properties  of  the  profile  can  be  determined.  The 
third  design  procedure  has  not  been  worked  out  in  detail,  so  that  no 
calculations  have  been  made  in  with  this  procedure. 

On  the  other  hand,  some  calculations  have  previously  been  made  with 
the  first  two  design  procedures.  However,  those  calculations  were  carried 
out  only  to  illustrate  the  efficacy  of  the  theory  [2,  3],  No  attempt  has 
yet  been  made  to  put  forward  specific  profile  designs  or  to  use  these 
two  procedures  to  explore  the  outcome  of  various  design  philosophies. 

The  purpose  of  the  present  study  is  to  remedy  this  deficiency.  In  order 
to  update  the  numerical  procedures  of  References  2 and  3 so  that  the 
computer  program  is  compatible  with  the  IBM  370  computer,  the  entire 
program,  as  described  in  Reference  3,  has  been  reviewed,  all  necessary 
corrections  and  modifications  have  been  made,  and  all  analytical  results 
of  Reference  2 have  been  checked  to  be  sure  that  the  calculated  results 
are  accurate.  This  work  has  been  carried  out  by  R.  F.  Davis  and  J. 
Fernandez.  They  have  documented  their  work  in  a companion  report  [8], 

It  should  be  emphasized  that  the  present  study  explores  only  the 
consequences  of  various  hydrodynamic  design  strategies  on  profile 
performance  and  geometry.  The  only  consideration  relating  to  profile 
structural  requirements  is  expressed  in  terms  of  the  cavity  thickness. 
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This  hydrodynamic  variable  is  treated  parametrically  in  order  to  show 
systems  and  structural  engineers  possible  hydrodynamic  consequences  of 
structural  or  other  systems  requirements.  The  important  design  inter- 
actions between  detailed  structural  and  hydrodynamic  requirements  are 
outside  the  scope  of  this  study. 

Flow  Geometry  and  Nomenclature: 

Before  turning  to  the  study  of  numerical  results  from  the  first  an 
second  design  procedures,  we  will  define  the  flow  geometry  and  the  maj  r 
parameters  which  are  basic  to  our  considerations  with  the  help  of  Fig  re  1. 
Although  this  schematic  illustrates  the  way  one  might  load  a profile  :or 
calculations  in  accordance  with  the  first  design  procedure,  the  esse  :ial 
geometric  aspects  of  the  hydrofoil,  the  cavity  and  the  free-stream  f dw 
apply  to  both  the  first  and  the  second  design  procedures.  In  the  c i e of 
the  second  design  procedure  the  region  of  secondary  perturbation  pr  ssure 
between  the  primary  peak  and  the  trailing  edge,  shown  in  Figure  2,  would 
be  deleted.  As  we  have  already  noted,  the  second  procedure  make  use  of 
an  angle  of  attack  increment  above  that  needed  for  shockless  er  .ry  to 
provide  for  the  prescribed  cavity  thickness  T.  In  this  case  ae  complete 
pressure  distribution  consists  of  two  parts  and  the  first  p?  t,  due  to 
the  nose  singularity  resulting  from  the  incremental  attack  ; ngle  above 
that  for  shockless  entry,  provides  automatically  that  the  ’ '©file  wetted 
surface  will  have  a non-zero  perturbation  pressure  everywh-  ’e,  except  at 
the  trailing  edge  where  it  must  be  zero.  Thus,  one  needs  > prescribe 

s 

only  the  primary  pressure  peak  of  the  second  part  in  order  to  determine 
the  shape  of  the  wetted  surface  when  the  second  design  procedure  is  used. 
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As  Figure  1 shows,  the  free-stream  velocity  and  static  pressure  are 
U and  P^,  respectively.  The  cavity  pressure  is  P^,  a constant,  and  if 
the  density  of  the  liquid  is  p the  cavitation  number  for  the  flow  is 


K 


- P. 
k 

1/2  p U5 


(1) 


If  the  static  pressure  at  any  point  in  the  flow  is  P(x,y),  we  shall 
define  the  perturbation  pressure  as 


p(x,y)  = 


P(x,y)  - Pfc 
1/2  p U2 


(2) 


with  the  understanding  that  at  points  on  the  wetted  surface  of  the 
hydrofoil  we  have  p(x,o)  = p(x),  where  now  0 <;  x 1 1.  On  the  upper 
and  lower  surfaces  of  the  cavity,  p(x,y)  = 0.  The  angle  of  attack  is 
a,  the  cavity  length  is  1 and  the  peak  pressure  location,  measured  from 
the  nose  of  the  foil  is  s.  Note  also  that  the  hydrofoil  has  a unit  chord 
and  that  the  cavity  thickness,  T,  is  defined  at  the  point  x=l  only. 

Figure  1 also  shows  the  positive  senses  of  the  lift,  cavity  drag  and 
the  pitching  moment  referred  to  the  nose  of  the  hydrofoil.  The  term 
cavity  drag  is  used  to  emphasize  that  skin  friction  on  the  profile 
wetted  surface  is  excluded  from  the  present  considerations.  If  the 
Reynolds  number  for  the  flow  is  known  this  added  drag  can  be  estimated 
separately.  It  is  much  smaller  than  the  cavity  drag.  The  three  hydro- 
dynamic  coefficients  assume  their  customary  definitions  for  two  dimen- 


sional flow: 
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provided  the  chord  length  is  unity.  In  this  study  the  upper  contour  of 

the  cavity  is  denoted  by  y^(x) , the  wetted  surface  contour  (or  camber 

function  by  T|(x);  and  the  strength  of  the  nose  singularity  corresponding 

to  an  attack  angle  above  shockless  entry  by  A.  The  relationship  between 

the  strength  of  the  point  drag  singularity  E,  and  the  incremental  cavity 

thickness  provided  by  the  complementary  function,  y (1),  is 

c 


E = y (1)//£T 

C 


where  & is  the  cavity  length  as  illustrated  in  Figure  1.  This  notation 
is  consistant  with  that  used  in  References  [1],  [2]  and  [3].  Additional 
information  on  this  and  relf.ted  theoretical  topics  are  given  in  those 
sources . 


FIRST  DESIGN  PROCEDURE-SUMMARY  OF  RESULTS 
Preliminary  Remarks: 

In  order  to  explore  design  trends  to  be  obtained  with  the  first 
design  procedure  a series  of  inverse  calculations  have  been  made  which 
use  the  pressure  distribution  of  Figure  2.  In  these  first  calculations 
the  location  of  the  peak  pressure  is  restricted  mostly  to  values  of  s between 


-17- 


June  30,  1975 
BRP:RFD:JF: jep 


.1  and  .5  chord.  For  the  present  this  limitation  will  not  impose  any 
essential  drawback  to  this  preliminary  investigation.  More  complete 
correlations,  corresponding  to  the  range  .lisi.9,  are  discussed  later 
in  this  report,  and  a tabulation  of  coefficients  for  the  entire  range  is 
given  in  Appendix  A.  The  range  of  CT  values  in  this  tabulation  has  been 
chosen  to  correspond  to  those  likely  to  be  encountered  in  the  design  of 
cavitated  foils  for  a high  speed  hydrofoil  boat. 

Figure  2 shows  a double  ellipse  pressure  distribution  having  a 
peak  pressure  magnitude  of  h and  a secondary  peak  of  Xh  where  0<X«1. 

The  secondary  ellipse  is  used  simply  to  insure  that  the  disturbance 
pressure  on  the  profile  will  exceed  zero  everywhere  on  the  wetted  surface 
of  the  profile  except  at  the  nose  and  at  the  trailing  edge.  The  ellip- 
tical distribution  has  the  form  p (x) Vl -x  near  x-1  in  accordance  with 
function-theoretical  requirements  appropriate  to  the  boundary  value  problem 
[2] . For  this  pressure  distribution  the  primary  semi-ellipse  always 
starts  at  the  point  x®0  for  all  values  s of  the  horizontal  semi-axis  for 
which  O^si.5  as  indicated  in  Figure  2.  In  those  cases  for  which  .5<s£l., 
the  semi-axis  of  the  primary  elliptical  distribution  is  (1-s)  and  this 
loading  is  anchored  to  the  trailing  edge  for  all  values  of  s corresponding 
to  this  rearward  loading  of  the  profile.  Thus,  the  primary  semi-ellipse 
occupies  the  interval  (2s-l)<x<l  and  the  quarter  ellipse  of  the  secondary 
loading  would  occupy  the  interval  0£x£(s2-l).  That  is,  if  we  reflect  the 
contours  of  Figure  2 in  the  vertical  axis  so  that  the  positions  of  nose 
and  trailing  edge  are  reversed  we  will  have  the  loading  configuration  for 
tail-loaded  profiles.  If  C^,  * and  X are  prescribed,  it  follows  at 
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once  that 


2C„ 


h - 


ir[s  + X(l-2s)l  * 


0<sil/2 


(7) 


When  l/2^s<l  one  has  the  case  of  aft  loading  on  the  foil  and 


2C. 


tt[1-s  + X(2s-1)]  * 


(8) 


Equation  (7)  determines  the  value  of  h for  the  various  lift  coefficients 
considered  in  the  initial  phases  of  our  study  of  nose-loaded 
profiles . 

Hydrodynamic  Performance : 

The  best  value  of  X to  be  used  for  foil  design  is  somewhat  uncertain 

because  the  designer  needs  only  enough  positive  pressure  in  the  secondary 

ellipse  to  provide  him  some  assurance  that  premature  wetted  surface 

cavitation  will  not  occur  within  some  range  of  attack  angles  less  than 

the  ideal  design  value  in  any  particular  case.  On  the  other  hand  if  X 

is  too  large  then  it  becomes  less  convenient  for  one  to  shift  the  center 

of  pressure  location  to  whatever  value  may  be  desired.  Clearly,  X»0  is 

the  best  choice  in  this  latter  respect.  In  order  to  examine  the  effect 

of  the  value  of  X a series  of  calculations  were  made  at  zero  cavitation 

number  and  fixed  cavity  thickness,  T*=0.15,  for  Xo,  .025  and  .05.  Values 

of  s were  selected  at  s“.3,  .4  and  .5.  Values  of  C.  were  chosen  to  be 

L 

at  » .08,  ,10,  .12,  .15,  The  important  parameter  for  performance 
considerations  is  the  lift-to-drag  ratio.  The  result  of  these  calcula- 
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tions  is  shown  in  Figure  3 in  which  L/D  versus  s curves  are  given  for 
X=0  and  X=.05.  The  curves  for  X=.025  lie  between  the  curves  for  X=0 
and  X=,05  and  they  were  omitted  for  the  sake  of  clarity.  For  future 
reference  this  illustration  also  shows  a curve  obtained  from  calculations 
made  with  the  second  design  procedure.  The  first  point  to  be  noticed 
from  this  illustration  is  that  a value  of  X as  large  as  .05  does  not 
affect  the  value  of  L/D  too  drastically  and  yet  it  should  provide  some 
margin  for  error  in  the  operating  attack  angle.  Therefore,  the  value 
X-.05  has  been  selected  as  the  standard  value  for  use  in  subsequent 
calculations  by  the  first  design  procedure. 

Another  trend  shown  by  Figure  3 is  the  increase  in  L/D  as  the  peak 
pressure  moves  toward  the  nose  for  all  C values.  This  is  an  important 
point  to  note  because  in  the  past,  designers  have  often  tried  to  place 
the  peak  pressure  as  near  to  the  trailing  edge  as  possible  in  order  to 
obtain  foils  of  low  drag.  This  design  philosophy  originated  with  Tulin 
[9]  and  it  was  further  reinforced  by  the  more  general  considerations 
of  Reference  [2] . However,  those  recommendations  ignore  the  fact  that 
the  upper  surface  of  the  cavity  must  clear  the  wetted  surface  of  the 
hydrofoil.  In  the  present  calculations  this  fact  is  explicitly  accounted 
for  in  the  design  procedure  and  it  results  in  a trend  with  respect  to 
the  location  of  the  peak  loading  which  is  directly  contradictory  to  the 
previous  result  — but  at  the  cost  of  less  favorable  lift-to-drag  ratios. 

A further  exploration  of  these  preliminary  trends  shows  the  effect 
of  cavitation  number  on  lift-to-drag  ratio  when  the  cavitation  number 
varies  from  0 to  0.2.  The  cavity  thickness,  T».15,  and  location  of  peak 
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pressure  at  s».3  were  held  fixed.  The  same  values  of  CT  were  selected 
as  before.  The  result  of  these  calculations  is  shown  in  Figure  4. 
Evidently  the  lift-to-drag  ratio  is  strongly  affected  by  cavitation 
number  with  the  most  favorable  values  being  found  at  K=0. 

Figure  5 shows  the  lift-to-drag  ratio  as  a function  of  peak  pressure 
location  for  various  lift  coefficients  at  a cavitation  number,  K“,2. 

Note  that  the  cavity  thickness  for  this  case  ia  T“.l  so  that  for  s*.3 
one  can  compare  this  illustration  with  Figure  4 and  see  that  increasing 
T at  K«=>.2  increases  the  cavity  drag.  This  trend  is  obviously  true  at 
all  K values.  Moreover,  we  see  from  these  curves  that  when  K".2  and 
T**.10  the  best  position  for  the  peak  pressure  varies  with  C . For 
example  at  Cj“.15  the  best  value  of  s appears  to  be  at  s«.42  and  that 
it  apparently  moves  toward  the  trailing  edge  as  decreases. 

In  Figure  6 the  lift-to-drag  ratio  has  been  plotted  against  peak 
pressure  location,  s,  at  a constant  Cj®.15  and  at  constant  T-.15  for 
various  cavitation  numbers  between  K**0  and  .2.  One  con  comparo  the 
curve  for  K*,2  with  the  curve  for  CL«.15  of  the  preceding  figure  to 
ser  again  how  increasing  the  cavity  thickness  from  T-.10  to  T-.15 
decreases  the  lift-to-drag  ratio  from  about  8.5  to  about  6.  For  the 
thicker  cavity  considered  in  Figure  6 it  appears  that  the  most  favorable 
location  for  the  peak  pressure  is  near  the  profile  nose. 

The  foregoing  results  from  the  first  design  procedure  have  revealed 


the  following  trends.  For  most  cases  one  will  obtain  sore  favorable 
lift-to-drag  ratios  if  he  concentrates  the  pressure  loading  as  near  to 
the  leading  edge  of  the  profile  as  he  can.  The  cavitation  number  has  a 
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large  effect  on  the  lift-to-drag  ratio.  The  larger  the  value  of  K 
becomes,  the  less  favorable  is  the  resulting  lift  to  drag  ratio.  This 
is  a direct  consequence  of  the  fact  that  one  requires  a fixed  value  of 
cavity  thickness,  T,  at  the  design  point  regardless  of  the  cavitation 
number.  But  as  the  cavitation  number  increases  the  cavity  length 
decreases  and  if  left  to  itself  the  cavity  thickness  would  tend  to 
decrease  too.  However,  our  design  method  requires  a constant  cavity 
thickness  at  the  trailing  edge.  Therefore,  it  is  necessary  that  we 
increase  the  strength  of  the  point  drag  function  to  meet  this  require- 
ment. Accordingly,  the  lift-to-drag  ratio  must  decrease  because  the 
cavity  drag  is  increased  while  we  hold  the  lift  constant.  This  situation 
should  be  true  for  all  locations  of  the  peak  pressure.  The  preceding 
examples  also  show  that  the  design  lift  coefficient  has  an  appreciable 
effect  on  the  lift-to-drag  ratio.  For  noso-loadod  profiles  the  lift-to- 
drag  ratio  increases  with  increasing  values  of  lift  coefficient.  This 
effect  seems  to  become  stronger  as  the  peak  pressure  location  is  moved 
toward  the  nose.  Now  it  Is  known  for  a delta  function  pressure  at  the 
trailing  edge  (2,  1]  that 


L ^ 8n(14K)  £(£-1) 
® (£-2)2 

where 


(9) 


vE-1 


4*(1+K)B 


(10) 


and 
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As  K-K),  and  equation  (9)  reduces  to  Tulin's  result  [9]  in 
which  L/D  = 8tt/C.  . The  important  thing  to  note  is  that  these  formulae 
indicate  that  the  lift-to-drag  ratio  should  decrease  as  increases 

which  is  opposite  to  the  trend  found  by  means  of  the  first  design 
procedure.  Thus,  if  structural  limitations  allow  it,  the  present  study 
suggests  that  the  designer  should  use  heavily  loaded  foils  in  order  to 
secure  the  most  favorable  values  of  L/D  at  the  design  condition.  The 
"delta  function  rule"  suggests  that  lightly  loaded  foils  are  preferred. 

We  have  already  noted  that  the  delta  function  rule  does  not  control  the 
position  of  the  cavity  upper  surface  with  respect  to  the  wetted  surface 
of  the  profile.  Evidently,  the  control  of  the  cavity  location  incor- 
porated into  the  present  design  procedure  changes  the  direction  of  the 
trend  of  L/D  with  respect  to  design  in  a most  significant  way.  However, 
we  will  sec  that  as  the  location  of  the  peak  pressure  moves  nearer  to 
the  trailing  edge  the  increase  of  L/D  with  increasing  becomes  less, 
although  this  trend  is  not  reversed  for  profiles  defined  by  the  first 
design  procedure.  Another  general  trend  resulting  from  the  foregoing 
calculations  is  tliat  as  the  cavitation  number  increases,  changes  in  the 
peak  pressure  location  have  less  effect  on  the  value  of  L/D.  Per  example. 
Figure  6 shows  tliat  the  lift-to-drag  ratio  of  nose  loaded  profiles  In 
the  range  .lfsi.5  is  hardly  changed  by  the  location  of  the  peak  pressure 
when  K*=  .2,  although  the  change  is  appreciable  when  K*0. 
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Hydrodynamic  Equivalence  and  Performance  Summary: 

All  of  the  preceding  analysis  of  trends  have  been  made  for  various 
peak  pressure  locations,  denoted  by  the  parameter  s of  the  double  ellipse 
pressure  distribution.  While  this  has  been  a convenient  procedure,  the 
foil  loading  quantity  of  fundamental  interest  is  x,  the  chordwise  distance 
from  the  hydrofoil  nose  to  the  center  of  pressure.  Values  of  x for  all 
cases  calculated  for  the  double  ellipse  pressure  distribution  are  tabulated 
in  Appendix  A.  Now,  cs-  has  considerable  freedom  in  the  specification 
of  the  pressure  disL.j.ibution  for  calculations  in  accordance  with  the 
first  design  p ccedure.  It  is  certainly  possible  to  prescribe  several 
distributions  which  will  have  the  same  value  of  x even  though  their 
shapes  may  differ  in  detail.  One  would  anticipate  that  as  long  as  all 
of  these  distributions  have  the  same  value  of  x,  CL,  T and  K.  at  the 
design  point,  the  remaining  hydrodynamic  performance  design  parameters 
should  differ  very  little  from  one  another.  In  order  to  test  this  con- 
jecture a series  of  calculations  have  been  made  with  the  first  design 
procedure  for  a modified  three-term  [2]  and  a reversed  three-term  pressure 
distribution  at  a C = .12,  as  shown  in  Figure  7.  The  reversed  three- 

Li 

term  distribution  has  x « ,3647,  corresponding  to  a nose-loaded  design. 

The  modified  three-term  distribution  has  x » .6685,  corresponding  to  a 
tail-loaded  design.  Hydrodynamic  performance  data  from  these  two  distri- 
butions are  compared  with  data  from  equivalent  double-ellipse  distributions 
in  Tables  I and  XI.  It  can  be  seen  that  differences  in  performance 
between  equivalent  profiles  are  very  small  indeed.  In  fact  we  believe 
that  part  of  the  differences  shown  may  be  due  to  the  fact  that  data  for 
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the  elliptical  distributions  were  obtained  from  those  of  Appendix  A 
by  linear  interpolation. 

The  foregoing  is  a useful  result  because  it  permits  one  to  summarize 
ax.L  of  the  hydrodynamic  trends  considered  above  for  the  double  ellipse 
family  of  pressure  distributions  with  the  understanding  that  this  summary 
will  also  be  valid  for  other  distributions  having  different  shapes.  All 
that  we  require  is  that  the  pressure  distributions  are  hydrodynamically 
equivalent  in  the  sense  that  the  value  of  x be  the  same  for  each  and  that 
the  same  design  point  (C^,  T and  K)  applies  to  all  of  them.  As  a result 
of  these  considerations  a series  of  "performance  maps"  have  been  prepared 
from  the  data  of  Appendix  A.  These  are  given  in  Figures  8,  9,  10,  11 
and  12,  which  show  plots  of  L/D  vs  with  contours  of  constant  and 
s.  Each  one  of  these  maps  shows  three  designs,  corresponding  to  three 
values  of  T,  and  each  map  is  valid  for  the  value  of  K noted  on  it. 
Evidently,  all  contours  are  equilateral  hyperbolas.  Note  that  the 
plot  for  T«,10  shows  intersecting  contours  for  certain  values  of  s. 

This  occurrence  is  a consequence  of  the  trends  shown  already  in  Figure 
5.  As  we  have  seen,  when  the  cavitation  number  is  high  enough  and  T is 
sufficiently  small,  the  best  location  for  the  peak  pressure  moves  from 
the  nose  toward  the  trailing  edge.  Evidently  the  nose  location  is  not 
always  the  best  location.  However,  this  is  the  only  case  found  so  far 
In  which  the  nose  location  is  not  the  best  choice.  The  chart  of  Figure 
13  is  a plot  of  vs  x,  with  contours  of  constant  and  s as  before. 
This  chart  provides  the  connection  between  x and  s so  that  if  the  design 
and  x are  known  the  resulting  value  of  s can  be  used  to  enter  the 
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charts  of  Figures  8 through  12  to  find  the  design  L/D,  etc.  for  distri- 
butions having  other  than  the  double  ellipse  shape.  Note  that  Figure 
13  applies  to  all  values  of  K and  T so  that  only  one  such  chart  is  needed. 
The  value  X,  the  secondary  quarter-ellipse  parameter,  for  these  charts 
is  .05. 

Profile  Geometry: 

The  foregoing  paragraphs  give  a fairly  comprehensive  overview  of 
the  hydrodynamic  performance  of  fully  cavitating  hydrofoil  sections 
designed  by  the  first  procedure.  It  is  now  necessary  to  explore  the 
geometric  properties  of  the  profiles  as  the  design  point  and  the  peak 
pressure  location  are  varied.  The  graphs  of  Figures  14  through  23  show 
to  scale  the  profile  wetted  surfaces  and  the  upper  cavity  contours  in 
accordance  with  the  nomenclature  defined  in  Figure  1.  In  these  figures 
the  chord  lines  of  all  profiles  are  on  the  x-axis  from  x=0  to  x»1.0. 

The  y-axis  shows  contour  ordinates  to  the  same  scale  as  the  abscissa  of 
each  point.  Note  that  only  that  portion  of  the  cavity  upper  surface 
which  lies  directly  above  the  wetted  surface  of  the  profile  is  plotted 
in  these  graphs.  The  trailing  edge  of  the  wetted  surface  is  located  at 
the  point  x®1.0  on  all  plots  and  the  lower  surface  of  the  cavity,  which 
would  extend  beyond  the  trailing  edge  Is  not  shown  at  all.  The  cavity 
uiickness,  T,  for  each  design  can  be  measured  directly  from  the  graphs 
ah  t.  vertical  distance  from  tLa  trailing  edge  of  the  profile  (1.0,  0) 
to  the  upper  surface  of  the  cavity  (1.0,  T).  The  space  between  these 
two  contours  from  the  profile  nose  at  (0,  0),  and  on  into  the  cavity 
beyond  the  trailing  edge  of  the  wetted  surface  if  need  be,  is  available 
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for  the  hydrofoil  structure.  We  have  not  drawn  in  any  foil  structure 
in  these  plots,  preferring  to  restrict  our  attention  to  the  hydrodynamics 
of  these  flows.  Of  course,  we  assume  that  there  would  be  some  clearance 
space  between  the  upper  surface  of  the  cavity  and  the  uppermost  parts  of 
the  profile  structure. 

The  profiles  of  Figures  14,  15,  16,  17  and  18  have  all  been  designed 
for  T=.15,  K=0,  C=.15.  Values  of  s vary  from  s*.2  to  s=.8.  (Note  that 

It 

s should  be  listed  as  an  input  quantity  in  these  figures.)  This  sequence 
of  profiles  shows  how  the  maximum  camber  of  the  profile  moves  toward  the 
trailing  edge  as  the  peak  pressure  point  moves  toward  the  trailing  edge. 
Moreover,  as  the  peak  pressure  point  moves  to  the  rear,  the  rounding  of 
the  profile  nose  increases,  corresponding  to  an  increase  of  the  point- 
drag  singularity  at  the  nose.  This  increase  is  required  in  order  to 
maintain  the  value  T=.15  as  s increases,  thereby  shifting  the  peak 
pressure  to  the  rear.  It  can  be  seen  that  the  cavity  drag  almost  doubles 
as  s increases  from  .2  to  .8.  In  the  following  graphs  we  will  explore 
the  effect  on  profile  geometry  when  s is  held  fixed  and  the  design  point 
is  varied. 

The  effect  of  varying  cavity  thickness,  T,  when  s=>.5,  K“0  and 
«.15  is  shown  by  Figures  17,  19  and  20  for  T=.10,  .15  and  .20  respec- 
tively. Notice  that  the  cavity  drag  nearly  triples  and  that  the  profile 
nose  becomes  more  rounded  as  T increases  over  this  range. 

Comparison  of  Figures  17  and  21  shows  the  effect  of  increasing  the 

cavitation  number  from  0.0  to  0.2  at  the  fixed  values  of  T“.15,  C. °.15 

L 

and  s».5.  We  see  that  the  shorter  cavity,  denoted  by  L in  the  figures, 
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corresponding  to  the  higher  value  of  cavitation  number  requires  an 
increase  of  the  point-drag  function  to  maintain  a thick  cavity  and  this 
results  in  greater  rounding  of  the  profile  nose. 

Figures  17  and  22  show  the  effect  of  reducing  from  .15  to  .08 
respectively.  Notice  that  as  C decreases  the  amount  of  undercamber 

Jj 

decreases  and  the  nose  curvature  seems  only  slightly  affected.  In  this 
comparison  the  fixed  parameters  are  T*=.15,  K-0.0  and  s=.5. 

Figure  23  illustrates  the  kind  of  a design  one  produces  when  he 
specifies  a thick  cavity,  fairly  high  C^,  forward  loading  and  high  cavi- 
tation number.  As  also  indicated  from  the  performance  map,  the  lift- 
to-drag  ratio  is  not  very  good.  As  a general  comment,  it  has  been 
found  in  this  series  of  calculations  that  profiles  having  a far-foward 
loading  or  a far -rearward  loading  such  as  s*.l  or  s=.9  require  a higher 
point  density  over  the  pressure  distribution  so  that  the  profile  wetted 
surface  will  be  accurately  calculated.  This  behavior  results  from  the 
sharp  pressure  peaks  which  require  fairly  dense  point  distribution  for 
an  adequate  number  of  points  in  the  steep  regions  of  the  distribution. 
Less  "peaky"  pressure  distributions  can  be  defined  with  fairly  good 
precision  by  50  equal  intervals  between  x=0  and  x=l.  The  contours 
shown  in  this  report  generally  result  from  this  spacing  and  close  ' 
inspection  of  Figures  14,  19  and  23  which  have  s-values  of  .2,  .5  and 
.8  respectively,  show  the  onset  of  the  behavior  noted  above  when  only 
50  points  are  used. 

Finally,  we  note  that  the  region  of  under camber  on  the  wetted 
surface  of  all  designs  illustrated  is  quite  dependent  on  the  position 
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of  the  peak  pressure.  For  example,  if  one  consic  rs  a nose-loaded 
profile  and  he  studies  the  region  between  x=0  and  x=i2s,  he  will  note 
that  the  undercamber  is  concentrated  in  this  re  ion.  Of  course,  the 
nose  rounding  required  by  the  point-drag  funct ion  modifies  this  behavior 
in  the  vicinity  of  the  nose.  However,  near  t.  i point  x=2s  the  curvature 
of  the  wetted  surface  reverses  and  beyond  thir  point  the  wetted  surface 
shows  no  undercamber.  In  the  case  of  tail-J jaded  profiles  the  situation 
is  the  same  except  that  the  influence  of  th  point-drag  function  is 
farther  away  and  one  can  now  make  the  same  .eneral  observations  if  he 
replaces  the  interval,  0£x£2s  by  2(l-2s)£  £1*  It  appears  therefore, 
that  while  the  hydrodynamic  performance  o.  the  profile  depends  primarily 
on  the  location  of  the  center  of  pressu?  j,  the  geometry  of  the  wetted 
surface  is  sensitive  to  the  peak  pressu.  a location  and  other  details  of 
the  pressure  distribution  shape. 

SECOND  DESIGN  PROCEDURE  - SUMMARY  OF  RESULTS 
Preliminary  Remarks: 

The  exploration  of  design  trenc  ■ derived  from  the  second  design 
procedure  start  with  the  modified  p assure  distribution  shown  in  Figure 
24.  These  pressure  distributions  re  obtained  from  those  defined  by 
Equations  (7)  and  (8)  simply  by  p ating  X«=0.  As  we  have  already 
explained,  the  secondary  pressuri  distribution  which  is  used  in  the 
first  design  procedure  to  insure  that  p(x)>0  except  at  the  leading  edge 
and  trailing  edges  of  the  profi  e is  not  needed  for  this  case.  In  the 
second  design  procedure,  this  equirement  is  automatically  satisfied  by 
the  added  pressure  distributi  a resulting  from  the  departure  of  the 
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attack  angle  from  the  ideal  angle. 

Hydrodynamic  Performance: 

As  before,  the  important  performance  parameter  for  design  purposes 

is  the  lift-to~drag  ratio.  A series  of  calculations  made  in  accordance 

with  the  second  design  procedure  has  given  the  results  plotted  in  Figure 

25.  This  figure  shows  the  effect  of  peak  pressure  location  on  lift-to- 

drag  ratio  at  K=0  for  a range  of  section  lift  coefficients  and  for  two 

cavity  thicknesses.  Note  the  similarities  between  the  trends  shown  in 

this  graph  and  those  obtained  from  analogous  calculations  by  the  first 

design  procedure  in  Figure  3.  Note  also  that  the  range  of  C and  T 

L 

values  in  Figure  25  is  somewhat  narrower  than  those  used  in  our  explor- 
ation of  the  first-design-procedure.  This  narrowed  scope  is  caused  by 

the  fact  that  C and  T cannot  always  be  specified  independently.  Reasons 
1* 

for  this  behavior  have  been  discussed  in  a more  quantitative  but  somewhat 
different  way  in  Reference  1.  For  the  present  discussion,  it  is  suffi- 
cient to  note  that  if  is  fixed  and  T is  also  prescribed,  then  situations 
can  occur  in  which  the  value  of  a required  to  meet  the  required  cavity 
thickness  will  cause  higher  lift  than  the  specified.  The  use  of  the 
point  drag  function  as  in  the  third  design  procedure  is  needed  to  over- 
come this  difficulty.  Thus,  Figure  25  shows  that  the  lowest  CT  value 
which  could  be  attained  when  T».l  is  about  .12  and  when  T=.15  the  lowest 
is  about  .15.  When  T=.20  it  was  found  that  the  lowest  value  of 
exceeds  .16, 

In  spite  of  the  limitations  discussed  above,  one  can  still  draw 
useful  conclusions  about  the  trends  produced  by  the  second  design 
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procedure  and  these  trends  can  also  be  compared  with  corresponding  results 
from  the  first  procedure.  Of  course  we  have  already  seen  in  the  case 
of  the  first  design  procedure  that  profile  performance  relates  primarily 
to  center  of  pressure  and  that  while  the  center  of  pressure  is  related 
to  peak  pressure  location,  s,  this  parameter  tends  to  have  an  important 
influence  on  wetted  surface  geometry.  Therefore  comparisons  of  perfor- 
mance trends  from  the  two  design  procedures  in  terms  of  s are  suggestive 
only.  It  remains  to  be  seen  how  they  compare  on  the  basis  of  center  of 
pressure  location. 

For  example,  let  us  compare  the  curve  for  C =.15  and  T=.15  from  Figure 
25  with  corresponding  curves  of  the  first  design  method  in  Figure  3 above. 
This  comparison  suggests  that  for  nose  loaded  foils  (s<.3)  the  first  design 
procedure  provides  profiles  having  more  favorable  lift-to-drag  ratios.  On 
the  otherhand,  when  the  peak  pressure  moves  toward  the  trailing  edge,  the 
second  design  procedure  is  better.  In  fact.  Figure  8,  the  first  design 
procedure  performance  map  for  K=0,  shows  that  L/D=7.09  when  s=.9,  T=.15 
and  Cj=.15.  This  is  only  about  60%  of  the  corresponding  L/D  value  at 
s=.9  in  Figure  25.  Thus,  when  location  of  the  peak  pressure  is  the  chief 
concern  and  a tail-loaded  profile  must  be  used  (s<.4)  the  second  design 
procedure  is  probably  the  better  design  procedure.  However,  if  securing 
the  most  favorable,  L/D  is  one's  primary  goal  then  it  would  seem  that 
nose  loading  is  to  be  preferred  and  the  first  design  procedure  should  be 
employed.  In  order  to  reinforce  this  point.  Figure  26  compares  L/D  for 
nose-loaded  profiles  for  s=.2,  with  T=.10  and  .15,  for  a range  of  values. 
The  curves  for  s=.3  are  included  in  Figure  25  simply  for  reference. 

If  for  some  reason  one  needs  to  design  a tail-loaded  profile 
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in  spite  of  the  L/D  penalties  noted  here,  the  elliptical  pressure 

distribution  may  not  be  the  best  choice.  For  sufficiently  high  values 

of  CT  and  s,  turbulent  separation  could  be  produced.  A means  of  avoiding 

this  occurrence  for  profiles  which  use  the  first  design  method  is  given 

in  Appendix  C.  We  note  in  passing  that  Appendix  C illustrates  how  one 

can  avoid  separation  but  still  have  the  steepest  possible  pressure  rise. 

Should  it  be  needed,  this  approach  aims  to  prescribe  appropriate  pressue 

distributions  without  requiring  extensive  turbulent  boundary  layer 

calculations  and  a turbulent  separation  model,  thus  avoiding  the  sort  of 

calculations  described  in  Reference  11,  for  example. 

The  foregoing  comparisons  summarize  the  chief  design  factors  to 

be  observed  in  connection  with  prescribed  peak  pressure  locations  for 

the  second  design  procedure.  Other  performance  trends  parallel  those 

of  the  first  procedure  at  least  qualitatively.  This  fact  is  illustrated 

within  the  parameter-range  limitations  of  the  second  procedure  by  the 

performance  map  of  this  method  of  profile  design.  But  first,  the  matter 

of  hydrodynamic  equivalence  in  the  sense  defined  above  for  the  first 
* 

foil  will  be  examined  for  the  present  case  of  the  second  design  procedure. 

Two-Term  Profiles,  Hydrodynamic  Equivalence  and  Performance  Maps: 

As  a preliminary  to  the  main  body  of  the  theory,  Reference  1 contains 

a discussion  of  hydrofoils  at  zero  cavitation  number  and  related  airfoils 

which  extends  the  ideas  originated  by  Tulin  [9].  As  an  example  of  optimum 

design,  a "two-term"  profile  is  found  which  has  the  most  favorable  lift- 

to-drag  ratio  and  which  has  a cavity  which  clears  the  wetted  surface. 

_ . 

p.  24. 
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The  cavity  clears  the  wetted  surface  because  the  pressure  distribution 
on  the  foil  contains  an  angle  of  attack  term  besides  the  two  terms  due 
to  camber.  Moreover,  there  is  no  complementary  function  in  the  solution 
so  that  the  example  treated  resembles  the  more  general  second  design 
procedure.  As  a matter  of' fact  we  can  put  the  "optimum"  results  from 
Example  3 of  Reference  1 in  a form  which  permits  comparison  with  the 
present  results.  In  order  to  facilitate  this  comparison  and  to  make  this 
document  more  nearly  self-contained,  Appendix  D gives  an  excerpt  form 
Reference  1 which  contains  these  preliminary  considerations.  For 
example,  we  find  from  Equations  (16)  and  (36)  of  Reference  1 that 


C a _ a 

M 160  a 


and  from  Equations  (16)  that 


°D*f  “2 


and  C.  « 


10 


a 


(12) 


(13) 


at  the  condition  for  most  favorable  lift-to-drag  ratio.  Since  xC.  « C. 

L J 

it  follows  for  this  foil  that  x = .2946.  Moreover,  at  a prescribed 
value  of  CT  we  can  write  the  lift-to-drag  ratio  as 


L 

D 


2Cl  U 25' 


(14) 


The  corresponding  result  for  a flat  plate  is  L/D  » tt/2Cl.  We  can  also 
use  Equations  (39)  and  (40)  of  Reference  1 to  write  the  cavity  thickness 
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at  the  trailing  edge  (x=l)  as 


T-l. 5764a  (15) 

which  permits  one  to  determine  CT , C and  C„  for  any  desired  cavity 

L 1)  M. 

thickness.  Table  III  shows  two-term  profile  performance  for  three 
values  of  cavity  thickness. 

From  this  example  one  sees  that  the  designer  has  more  freedom  if 
he  uses  the  more  general  second  design  procedure  because  the  "optimum" 
condition  selected  in  Reference  1 results  in  a fixed  portion  of  the  total 
lift  which  is  carried  by  that  part  of  the  pressure  distribution  due  to 
camber.  For  example,  in  the  second  design  procedure  we  express  that 
fraction  of  CL  which  is  carried  by  the  prescribed  pressure  distribution 
as 


V “ m CL  * (16) 

Then  that  part  of  carried  by  the  extra  angle  of  attack,  or  the  effect 
of  a in  Example  3 of  Reference  1,  is  (l-m)CT . We  can  make  use  of  Equa- 
tions  (36)  and  (27)  from  Reference  1 to  determine  that 

m - 6/7  » (17) 

which  holds  for  all  values  of  a.  That  is  to  say,  m is  fixed  once  and 
for  all  and  this  result  is  consistant  with  the  fact  that  x is  also 
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independent  of  a for  this  particular  example  of  an  "optimum"  profile. 
This  is  not  the  case  for  the  second  design  procedure  as  can  be  seen 
from  the  data  of  Appendix  B.  In  the  second  design  procedure,  we 
specify  C^,  T and  K as  well  as  the  shape  of  the  pressure  distribution 
due  to  camber.  As  a result  the  parameter  m shows  considerable  variation 
within  its  permissible  range,  0£m£l. 


In  spite  of  the  differences  just  discussed,  it  seems  useful  to 


compare  the  performance  of  a two-term  profile  with  an  equivalent  profile 


obtained  from  the  second  design  procedure.  To  do  this  we  compare  the 


results  of  Table  III  with  and  T values  in  Appendix  B for  K-Q.  We 
see  that  if  T is  greater  than  .10,  the  values  of  CL  in  Table  III  are 
outside  the  range  of  cases  studied  in  Appendix  B.  However,  when  T>.10 


“ .1395.  But  =>  .14  is  one  of  the  lift 
coefficients  tabulated  in  Appendix  B.  These  values  differ  by  1.1Z 
and  if  one  considers  the  effect  of  this  variation  on  the  data  of 

Appendix  B for  s«.3,  he  will  conclude  that  L/D  would  suffer  a corresponding 
variation  of  about  0.5X.  Therefore,  it  seems  reasonable  to  interpolate 
for  that  value  of  s corresponding  to  x » .2946.  This  corresponding  value 
is  8®. 2834.  Other  quantities  relating  to  both  two-term  and  elliptical 
pressure  distributions  are  compared  in  Table  IV.  From  this  table  one 
can  conclude  for  K»0,  at  least,  that  shape  variations  among  hydrodynomi- 
cally  equivalent  pressure  distributions  will  have  little  effect  on  the 
hydrodynamic  performance  of  the  resulting  profiles.  Note  that  tho 
value  of  m for  the  two  cases  shows  a significant  difference. 

Of  course,  the  foregoing  results  apply  strictly  to  tho  second  design 
procedure  at  zero  cavitation  number.  These  findings  need  to  be  extended 


the  two-term  foil  has 
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to  other  values  of  the  cavitation  number  such  as  K-.05  and  K=,10  and  to 
other  cavity  thicknesses,  such  as  T=,10  and  T=.15.  A convenient  choice 


for  the  lift  coefficient  is  «*  .16  because  Appendix  B shows  a good 


range  of  x at  both  values  of  K.  For  the  comparison  calculations,  we 
will  use  the  modified  three-term  and  reversed  three-term  pressure 
distributions  of  Figure  7 except  that  the  ordinates  of  that  illustration 
must  be  scaled  in  the  ratio  16/12  in  order  to  give  the  design  in  the 
event  that  m*l.  The  results  of  all  calculations  are  given  in  Tables 
V and  VI.  As  in  previous  comparisons,  the  performance  of  the  equivalent 
elliptic  pressure  distribution  is  obtained  by  linear  interpolation  from 
values  in  Appendix  B. 

Tables  V and  VI  show  again  that  the  performance  of  hydrodynamically 
equivalent  profiles  is  insensitive  to  detailed  variations  of  pressure 
distribution  shape.  An  interesting  aspect  shown  by  these  tabulations 
and  Table  IV  also  is  that  the  m-values  for  the  ollipticol  distribution 
arc  always  somewhat  smaller  titan  those  for  the  three-term  or  the  two-term 
profiles.  As  was  found  in  the  case  of  the  first  design  procedure,  it 
also  appears  in  the  case  of  the  second  design  procedure  that  profile 
performance  estimates  can  be  based  on  data  from  the  elliptic  pressure 
distribution  as  long  as  the  pressure  distributions  are  hydrodynamically 
equivalent.  The  practical  problem  for  the  application  of  this  rule  to 
the  second  design  procedure,  is  that  one  new  needs  to  carry  out  the 
design  calculation  for  the  new  profile  in  order  to  establish  the 
equivalence  between  it  and  one  having  an  elliptic  pressure  distribution. 
This  situation  is  caused  by  the  fact  that  now  the  pressure  distribution 
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has  two  parts  and  we  can  specify  only  that  part  which  affects  the  camber 
of  the  profile  whereas  the  center  of  pressure  coordinate  depends  upon 
the  relative  magnitude  of  both  parts.  This  fact  prevents  us  from  using 
a summary  of  results,  such  as  that  presanted  for  the  first  design  pro- 
cedure in  the  charts  of  Figures  8 through  13,  as  rather  precise  guides 
for  profile  design  performance. 

On  the  other  hand,  performance  maps  which  summarize  the  properties 

of  profiles  having  elliptic  pressure  distributions  do  provide  a convenient 

overview  of  trends  to  be  expected  from  the  second  design  procedure  and  they 

can  give  the  designer  qualitative  guidance.  Figures  27,  28  show  the  forces 

and  Figures  29,  30  and  21  show  the  pitching  moment.  These  figures  illustrate 

the  results  tabulated  in  Appendix  B,  except  that  moment  data  for  T=.15  have 

not  been  plotted  in  order  to  simplify  these  illustrations.  As  discussed 

previously,  the  range  of  parameters  for  the  second  design  procedure  is  not 

as  extensive  as  it  is  for  the  first  design  procedure  and  so  there  are  fewer 

charts  of  L/D  versus  in  this  case.  However,  there  are  more  charts  of 

Cy  versus  Sc  because  uov  these  contours  depend  upon  K and  T. 

Center  of  Pressure  Location  and  Comparison  of  First  and  Second  Design 
Procedures i 

All  of  the  preceding  comparisons  have  used,  s,  the  position  of  the 
peak  prescribed  pressure  cm  the  wetted  surface  as  the  primary  foil  dosign 
par&mtcr.  Tills  lias  been  a convenient  choice  because  it  can  be  prescribed 
quite  independently  of  the  final  performance  parameters;  and,  as  we  have 
seen  in  the  case  of  foils  derived  from  the  first  design  procedure,  it  is 
Important  with  respect  to  the  wetted  surface  geometry.  However,  the 
hydrodynamic  performance  of  profiles  designed  by  both  procedures  lias 
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beer  found  to  depend  on  x,  the  location  of  the  center  of  pressure.  This 
parameter  is  certainly  related  to  the  value  of  s,  but  for  the  second 
design  procedure  the  added  pressure  distribution  caused  by  the  increase 
in  attack  angle  above  the  ideal  angle  makes  this  relationship  somewhat 
indirect.  It  remains,  therefore,  to  see  how  the  idea  of  hydrodynamic 
equivalence  applies  to  comparisons  between  the  first  and  second  design 
procedures. 

Figures  32,  33  and  34  present  plots  of  L/D  against  x for  profiles 
obtained  from  the  two  design  procedures.  The  points  on  the  curves  show 
various  s-values  given  in  Appendices  A and  B.  In  Figure  32,  those 
points  nearest  the  L/D  axis  on  the  curves  for  the  first  design  procedure 
correspond  to  s=.l  and  on  the  curves  for  the  second  design  procedure, 
the  points  plotted  nearest  to  the  L/D  axis  are  for  sa.2  when  CL=.12  and 
s=.3  when  C^-,15.  This  situation  is  also  true  for  Figures  33  and  34, 
except  that  now  the  C =.15  curves  for  the  second  design  procedure  start 
at  sa.2  also.  These  figures  show  that  the  use  of  x for  comparing  first- 
and  second-foil  performance  only  exhibit  rather  small  L/D  differences 
over  the  range  of  x for  which  they  are  comparable.  Of  course  these 
differences  are  much  larger  than  those  displayed  in  Tables  I through  VI 
above.  But  those  data  compare  foil  designs  within  a given  procedure  so 
that  relatively  minor  variations  of  pressure  distribution  shape  are 
considered.  When  we  compare  a profile  from  the  first  design  procedure 
with  an  equivalent  profile  from  the  second  design  procedure,  the  pressure 
distributions  for  the  two  cases  differ  significantly.  Therefore,  we 
expect  to  see  more  of  a difference  for  the  present  comparisons  than  we 
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found  above.  The  fact  that  these  differences  are  as  small  as  they  are 
is  noteworthy.  Finally,  we  note  that  if  one  accepts  the  differences  shown 
as  being  significant,  we  can  state  that  theoretically  the  first  design 
procedure  generally  leads  to  profiles  having  the  lower  drag.  The  practi- 
cal importance  of  this  slight  advantage  is  questionable. 

The  limited  range  of  center  of  pressure  travel  of  the  second  design 
procedure  compared  to  the  first  design  procedure  is  another  feature 
which  deserves  comment.  For  example,  when  K=0,  it  is  known  that  x = 5/16 
for  the  flat-plate  hydrofoil  at  all  angles  of  attack.  But  the  added 
pressure  distribution  due  to  attack  angle  found  in  the  second  design 
procedure  is  the  flat-plate  distribution.  Therefore,  its  presence 
inhibits  the  degree  by  which  the  designer  can  vary  x by  changing  the 
value  of  s.  Thus,  we  see  that  all  of  the  x values  for  the  second  design 
procedure  are  less  than  .5  even  though  s values  as  large  as  .9  are 
included  in  Figures  32,  33  and  34.  We  also  see  that  the  forward  travel 
of  x available  to  the  designer  is  similarly  limited  although  this 
restriction  is  not  os  severe  as  it  is  for  rearward  travel.  This  behavior 
explains  the  apparent  superiority  of  the  second  procedure  compared  to 
the  first  procedure  for  large  values  of  s noted  in  Figure  3. 

Geometric  Properties; 

Discussion  of  performance  trends  for  the  first  and  second  procedures 
are  now  complete.  The  remaining  aspect  of  the  second  design  procedure 
requiring  discussion  is  the  profile  geometry  obtained  by  this  method. 

The  necessary  results  are  contniued  in  Figures  35,  36,  37  and  38  which 
show  wetted  surface  and  upper  cavity  contours  for  four  representative 
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flows  at  zero  cavitation  number.  These  plots  conform  to  the  same  rules 

given  previously  for  the  first  design  procedure.  Figures  35  and  36  show 

how  the  undercamber  of  the  vetted  surface  is  increased  as  the  design  CT 

L 

increases  from  .12  to  .15.  This  trend  was  also  observed  in  the  case  of 
foils  resulting  from  the  first  design  procedure.  Figures  36  and  37  show 
how  the  region  of  undercamber  moves  toward  the  profile  trailing  edge  as 
the  peak  pressure  location  is  changed  from  s=.3  to  s=.7.  This  trend 
has  also  been  observed  for  contours  derived  from  the  first  design  method. 
Finally,  a comparison  of  Figures  37  and  38  illustrates  how  the  wetted 
surface  changes  for  s=.3  when  increases  from  .15  to  .16  and  at  the 

same  time  T increases  from  .10  to  .15.  It  does  not  seem  necessary  to 
show  the  effect  on  wetted  surface  shape  of  increasing  K above  zero  while 
other  parameters  are  fixed  because  this  trend  is  essentially  the  same 
as  that  illustrated  for  the  first  design  method.  In  fact,  the  essential 
difference  between  the  two  procedures  as  far  as  foil  geometry  is  con- 
cerned is  that  the  leading  edge  of  profiles  obtained  from  the  second 
design  procedure  will  be  sharp  instead  of  rounded.  This  can  cause  some 
problems  structurally,  particularly  for  a situation  as  illustrated  in 
Figure  36,  In  fact  if  one  selects  s«.l  or  .2  and  does  not  require  a 
large  enough  value  of  T the  cavity  surface  will  intersect  the  wetted 
surface  contour.  This  occurrence  is  not  as  likely  with  the  first 
design  method.  However,  In  either  case  one  will  not  always  produce 
acceptable  designs  if  the  peak  pressure  is  moved  too  far  forward.  More- 
over, previous  remarks  concerning  computational  accuracy  needed  to  define 
the  vetted  surface  for  peak  pressure  locations  very  close  to  the  leading 
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edge  or  very  close  to  the  trailing  edge  also  apply  to  the  second  design 
procedure. 

OFF-DESIGN  CALCULATIONS 

Suppose  the  wetted  stir face  and  upper  cavity  contours  have  been 
determined  from  inverse  calculations  by  means  of  the  first  or  second  design 
procedure.  The  next  step  is  to  examine  the  off-design  performance  and 
critical  flow  geometries  of  the  candidate  hydrofoil  section.  In  this 
study,  we  shall  carry  out  such  calculations  by  means  of  the  linearized 
direct  theory  of  cavity  flows  [10].  Of  course,  we  know  that  linearized 
theory  lacks  the  precision  of  the  nonlinear  theory.  Nevertheless,  since 
the  present  inverse  theory  is  linearized  it  is  consistant  to  use  the 
linearized  direct  theory  with  it.  In  this  way,  one  obtains  a consistant 
set  of  preliminary  results  for  both  design  and  off-design  flow  states. 

Once  a few  promising  candidates  for  a particular  application  have  been 
found  from  the  linearized  theory,  one  can  complete  the  design  process 
by  recourse  to  nonlinearized  cavity-flow  calculations. 

Linearized  calculations  for  hydrofoil-cavity  interference  can  be 
carried  out  if  one  uses  the  direct  theory  of  Reference  8 to  determine 
the  upper  cavity  contour  at  any  attack  angle  a.  In  these  calculations, 
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a = a$  + W(a) 


(18) 


where  a = A-l  and  3 is  defined  by  Equation  (11) . The  precise  form  of 
all  terms  in  this  equation  is  given  in  Reference  2.  This  reference  also 
shows  that  the  upper  cavity  contour  at  any  off-design  condition  can  be 
expressed  as 


y(x)  = 1 (x,a)  + S2  (x,a) 


(19) 


where  the  functions  and  are  obtained  from  y(x)  in  Reference  2 by 
grouping  its  terms  so  that  $ is  a coefficient  of  as  shown. 

Suppose  next  that  one  selects  a contour  t(x)  which  lies  within  the  area 
bounded  by  the  wetted  surface  of  the  profile  and  the  upper  cavity  contour 
at  the  design  flow  state.  The  contour  t(x)  is  the  upper  surface  of  the 
foil  structure  will  simply  be  given  by  table  of  ordinates  over  the  interval 
Oixll.  Then  a critical  off-design  flow  state  will  exist  when  the  cavity 
contour  y(x)  is  tangent  to  t(x).  This  flow  state  can  be  found  if  one  puts 
y(x)“t(x)  in  the  off-design  cavity  contour  equation  above  and  solves  for  3: 


t (x)  - S2  (x,a) 
~z1  (x,a)  “ 


(20) 


Then  if  one  chooses  a value  of  a,  he  can  calculate  a value  of  3 for 
each  x in  (o,x).  He  can  select  the  largest  value  of  3 from  among  these 
values.  Then,  he  can  use  this  same  value  of  a and  this  greatest  value 
of  3 in  the  closure  equation  to  find  the  corresponding  a.  One  can  then 
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calculate  3 from  Equation  11.  The  three  values  (a,  K and  a)  define  an 
off-design  flow  state  at  which  the  cavity  is  just  tangent  to  the  upper 
(non-wetted)  surface  of  the  hydrofoil.  This  calculation  can  be  repeated 
for  a sequence  of  a-values  which  divides  an  a-K  plane  into  two  regions, 
in  the  region  above  this  curve  the  flow  will  be  at  a permissible  state. 
Points  on  or  below  this  curve  correspond  to  states  having  interference 
between  the  foil  and  the  cavity.  The  actual  equations  employed  and  the 
numberical  methods  used  in  these  calculations  are  given  in  Reference  3. 

An  example  of  results  obtained  from  the  calculations  described  above 
has  already  been  given  in  Reference  2.  Nevertheless  it  is  useful  to 
give  further  examples  which  correspond  to  flow  conditions  typical  of  this 
study.  Accordingly,  we  have  selected  the  parameters  s=.3,  CL=.12,  K=0 
and  T=.10  for  the  design  flow  state  for  two  profiles:  one  from  the  first 

and  one  from  the  second  design  procedure.  The  results  of  these  calcula- 
tions are  plotted  in  Figures  39  and  40  in  an  a-K  plane.  The  design  attack 
angle  for  these  profiles  is  plotted  on  the  o'-axis.  Points  in  this  plane 
above  the  curves  are  permissible  flow  conditions.  Those  which  are  on 
or  below  the  curves  are  not  permissible  flow  states.  In  these  figures 
"alpha  critical"  refers  to  the  locus  of  points  defining  the  conditions 
of  cavity  tangency.  In  these  calculations  the  upper  surface  of  the  foil 
has  been  derived  from  the  upper  cavity  surface  in  each  case  by  putting 

t(x)  = 6y(x)  , O<0<1  . (21) 

In  both  examples  we  have  put  6«.9  and  .8.  Thus,  the  hydrofoil  sections 
have  a wedge-like  shape. 
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To  continue  these  examples  we  have  also  included  the  general  off- 
design  calculation  method  of  Reference  8.  The  two  wetted-surface  slope 
functions  T)'(x)  noted  above  are  also  inputs  to  these  calculations  which 
are  carried  out  for  a range  of  K and  a values.  They  permit  us  to 
calculate  k,  CT , CL,  C..  and  L/D  for  the  prescribes  values  of  £ and  a. 
Results  for  the  two  profile  designs  are  shown  in  the  curves  of  Figures 
41  to  48.  It  will  be  noticed  that  we  have  also  superimposed  the  locus 
of  interference  conditions  from  Figures  39  and  40  when  8=.9  in  Figures 
41  and  45  to  show  how  interference  can  cut  off  the  available  range  of 
off-design  performance.  Aside  from  the  interference  data,  general  trends 
shown  in  these  illustrations  follow  those  which  have  been  noted  in  many 
previous  investigations  of  the  direct-theory  of  cavity  flows.  No 
further  discussion  of  those  found  here  seems  necessary. 

CONCLUSIONS 

Perhaps  the  most  important  general  conclusion  to  emerge  from  this 
study  is  that  the  most  favorable  lift-to-drag  ratios  are  most  often 
found  when  the  profile  is-  designed  to  have  its  center  of  pressure  placed 
as  near  to  the  leading  ed'  e of  the  profile'  as  possible.  Only  one 
exceptional  combination  of  T=,10  and  K=.  £ has  been  found  so  far.  This 
result  holds  for  profiles  designed  with  either  the  first  or  the  second 
design  procedure.  Of  course,  this  finding  must  be  applied  with  geometric 
and  other  modifying  factors  kept  in  mind.  Consequently,  this  result  will 
seldom,  if  ever,  be  applied  to  the  fullest  extent. 

For  foils  resulting  from  either  design  procedure,  it  is  found  that 
at  the  design  point  (prescribed  values  of  lift  coefficient,  cavitation 
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number,  and  cavity  thickness)  and  with  a pressure  distribution  fixed  in 


shape  but  having  variable  amplitude,  the  lift-to-drag  ratio  increases 


with  increasing  lift  coefficient.  That  is,  within  structural  limita- 


tions, the  profile  should  be  designed  to  be  as  heavily  loaded  as  possible. 


The  cavity  thickness  at  the  design  point  has  a powerful  effect  on 


the  lift-to-drag  ratio.  For  either  design  procedure,  as  the  the  thick- 


ness increases  the  cavity  drag  increases. 


The  cavity  thickness  at  the  design  point  also  influences  the  lift- 


to-drag  ratio  for  a given  design.  Increasing  the  cavitation  number 


decreases  the  lift-to-drag  ratio  for  foils  obtained  from  both  design 


procedures  for  most  attack  angles  of  practical  interest. 


Another  result  found  in  this  study  is  that  the  location  of  the 


center  of  pressure  is  the  most  important  property  of  the  pressure 


distribution  as  far  as  hydrodynamic  performance  is  concerned.  As  long 


as  the  center  of  pressure  is  held  fixed,  variations  in  the  shape  of 


the  pressure  distribution  have  almost  no  effect  on  the  hydrodynamic 


forces  on  the  profile.  Of  course,  fixing  the  center  of  pressure  does 


not  allow  for  very  large  variations  in  pressure  distribution  shape 


because  this  result  was  found  to  be  true  when  the  design  point  is  also 


fixed.  When  the  design  point  is  fixed  and  two  profiles  have  the  same 


center  of  pressure  they  are  said  to  be  hydrodynamically  equivalent. 


This  equivalence  is  particularly  precise  when  the  two  profiles  result 


from  the  same  design  procedure  but  have  differing  pressure  distribution 


shapes. 


Hydrodynamic  equivalence  has  also  been  found  to  be  a useful  concept 


' ' ' py ''''V  ' 
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when  one  profile  designed  by  the  first  procedure  is  compared  with  a foil 
designed  by  the  second  procedure.  If  both  have  the  same  design  point 
and  center  of  pressure  the  hydrodynamic  forces  are  nearly  the  same. 

While  this  equivalence  is  less  exact  than  it  is  for  profiles  obtained 
by  the  same  design  procedure  it  appears  that  agreement  of  the  forces  is 
close  enough  for  many  practical  purposes.  Of  course  fairly  small 
differences  were  found  in  the  forces  and  if  one  wishes  to  attach 
importance  to  them  he  could  say  that  the  first  design  procedure  leads 
to  slightly  more  efficient  profiles  than  the  second  procedure.  We  do 
not  believe  this  to  be  a very  important  point. 

Finally  it  was  found  that  the  shape  of  the  wetted  surface  is  most 
sensitive  to  the  location  of  the  peak  pressure.  Although  the  peak 
pressure  location  is  related  to  the  center  of  pressure  location,  its 
influence  on  the  undercamber  of  the  wetted  surface  appears  to  be  most 
direct  and  its  influence  on  profile  forces  seems  somewhat  less  direct 
than  does  the  center  of  pressure  location. 

The  limited  scope  of  this  study  suggests  that  although  we  have 

•k. 

established  desirable  design  goals  from  a purely  hydrodynamic,  viewpoint, 
structural  considerations  may  not  always  permit  their  fullest  exploitation. 
We  have  seen  that  enhancement  of  lift-to-drag  ratios  is  accomplished  by 
a reduction  in  the  cavity  cross-sectional  area  available  for  the 
profile  structure.  This  is  especially  true  near  the  leading  edge  where 
chordwise  and  spanwise  bending  stresses  in  a three  dimensional  hydrofoil 
may  become  too  high  to  meet  strength  requirements.  From  this  viewpoint 
it  might  be  better  to  compre  lift-to-drag  ratios  against  an  appropriate 
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strength  criterion  instead  of  using  the  parameters  of  cavity  thickness 
and  lift  coefficient.  Therefore,  the  present  study  is  not  a complete 
guide  for  the  use  of  inverse  theory  in  hydrofoil  design.  However,  we 
hope  that  the  hydrodynamic  trends  are  well  enough  established  to  make 
the  union  of  hydrodynamic  and  structural  design  considerations  less 
difficult. 
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Table  I - Performance  Comparison  for  Nose  Loading,  First  Design  Procedure 


K - .05 

1 

C 

1 D 

a0 

L/D 

.00704 

3.2808 

17.0557 

.00702 

3.2266 

• 17*.  0944 

.01338 

4.7132 

8.9672 

.01336 

4.6590 

8.9841 

.02174 

6.1456 

5.5201 

.02171 

6.0914 

5.5285 

Notes:  x = .3647 


L/D  C =.12,  C =-.04376 

L M 


Table  II  - Performance  Comparison  for  Tail  Loading,  First  Design  Procedure 


n 

K « .05 

K = .10 

Notes:  x = .6685 

D 

°D 

a0 

L/D 

CD 

a° 

L/D 

Cl=.12,  (^=-.08022 

.10 

.00914 

2.4476 

13.1299 

.01063 

2.4666 

11.2887 

3-TERM 

.00909 

2.4126 

13.2019 

.01060 

2.4331 

11.3224 

DOUBLE  ELLIPSE 

.15 

.01632 

3.8800 

7.3529 

.01838 

3.8990 

6.5291 

3-TERM 

.01626 

3.8500 

7.3842 

.01833 

3.8655 

6.5481 

DOUBLE  ELLIPSE 

.20 

.02551 

5.3124 

4.707 

.02812 

5.3313 

4.2675 

3-TERM 

.02542 

5.2774 

4.721 

.02805 

5.2979 

4.2788 

DOUBLE  ELLIPSE 
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Table  III  - "Two-Term"  Profile  Performance  at  Zero  Cavitation  Number 


T 

a0 

rv 

c» 

L/D 

C 

M 

.10 

3.6346 

.1395 

.00632 

22.07 

-.04110 

.15 

5.4519 

.2092 

.01422 

14.71 

-.06166 

.20 

7.2692 

.2790 

.02528 

11.04 

-.08221 

Table  IV  - Performance  Comparison  of  Equivalent  Twoj-Term  and  Elliptical 
Pressure  Dis  trlfcutions : T° . 10 „ K«0  and  x« . 2946 


.1395  . 

.1400  . 


% 

L/D 

a* 

PRESSURE  DISTRIBUTION 

0063 

22.07 

3.635 

TWO-TERM,  ra*.8571 

0064 

21,87 

3.602 

ELLIPTICAL,  a».6529 
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Table  V - Performance  Comparison  for  Aft-Loaded  Profiles,  Second  Design 
Precedure,  Design  C=.16 


-.0690  .4312 


-.0522  .3263 


.10  .10  -.0586  .3675 


CD 

L/D 

a0 

m 

PROFILE 

.0085 

18.73 

3.254 

.3331 

ELLIPTICAL 

.0086 

18.59 

3.255 

.3407 

3-TERM,  REVERSED 

.0142 

11.23 

5.126 

.3630 

ELLIPTICAL 

.0143 

11.22 

5.126 

.3709 

3-TERM,  REVERSED 

.0096 

16.61 

3.487 

.1398 

ELLIPTICAL 

.0096 

16.58 

3.488 

.1434 

3-TERM,  REVERSED 

S)  1 
U ! 

L/D 

a0 

m j 

PROFILE 

-.0540  .3375 


-.0506  .3162 


.10  ,10  -.0524  .3275 


,0076 

.0077 


.0141  11.33  5.157 
.0141  11.32  5.159 


.0093  17.15  3.607 
.0093  17.12  3.611 


0587  ELLIPTICAL 


.2247  I ELLIPTICAL 


2390  ' 3-TERM 
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Figure  1 - Schematic  Diagram  of  Two-Dimensional  Flow  Geometry 
and  Prescribed  Pressure  Distribution  for  the  First 
Hydrofoil  Design  Process 


Figure  2 - Double-Ellipse  Pressure  Distribution  for  Use  with 
the  First  Design  Procedure 


PERTURBATION  PRESSURE.  p(x) 
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Figure  7 - Three-Term  and  Reversed  Three-Term 

Pressure  Distributions  for  Comparison 
of  Results  with  Results  from  Elliptic 
Distributions 
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FIGURE  13.  PITCHING  MOMENT  PERFORMANCE  MAP 
FIRST  DESIGN  PROCEDURE 
DOUBLE  ELLIPSE  PRESSURE 
DISTRIBUTIONS,  X=  .05 
VALID  FOR  ALL  VALUES  OF  T AND  K 


Figure  14  — Cavity  and  Wetted  Surface  Geometry 

Double  Ellipse  Pressure:  Distribution,  X = .05,  First  Design  Procedure 


Figure  15  - Cavity  and  Wetted  Surface  Geometry 

Double  Ellipse  Pressure  Distribution,  X *=  .05,  First  Design  Procedure 


Figure  16  - Cavity  and  Wetted  Surface  Geometry 

Double  Ellipse  Pressure  Distribution,  X * .05,  First  Design  Procedure 


jure  17  - Cavity  and  Wetted  Surface  Geometry 

Double  Ellipse  Pressure  Distribution,  A = .05,  First  Design  Procedure 


'avity  nad  Wetted  Surface  Geometry 

Double  Ellipse  Pressure  Distribution,  X « .05,  First  Design  Procedure 


Figure  19  - Cavity  and  Wetted  Surface  Geometry 

Double  Ellipse  Pressure  Distribution,  X - .05,  First  Design  Procedure 


Figure  20  - Cavity  and  Wetted  Surface  Geometry 

Double  Ellipse  Pressure  Distribution,  X «*  .05,  First  Design  Procedure 


Figure  21  - Cavity  and  Wetted  Surface  Geometry 

Double  Ellipse  Pressure  Distribution,  X ■ .05,  First  Design  Procedure 


Figure  22  - Cavity  and  Wetted  Surface  Geometry 

Double  Ellipse  Pressure  Distribution,  A ■*  .05,  First  Design  Procedure 


Figure  23  - Cavity  and  Wetted  Surface  Geometry 

Double  Ellipse  Pressure  Distribution,  X **  .05,  First  Design  Procedure 
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Figure  26  - Further  Comparisons  of  Results  from 

Fir&t  and  Second  Design  Procedures  for 
Nose-leaded  Hydrofoils 
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FIGURE  29.  SECOND  DESIGN  PROCEDURE , P ITCHING 
MOMENT  PERFORMANCE  MAP 

CAVITATION  NUMBER,  K«0 

TsCAVITY  THICKNESS  AT  TRAILING 
EDGE  OF  PROFILE 

S*  DISTANCE  OF  PEAK  PRESSURE  FROM 
NOSE  OF  UNIT-CHORD  HYDROFOIL 
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FIGURE  30.  SECONO  DESIGN  PROCEDURE,  PITCHING 
MOMENT  PERFORMANCE  MAP 

CAVITATION  NUMBER,  K=. 05 

T= CAVITY  THICKNESS  AT  TRAILING 
EDGE  OF  PROFILE 

. 8 'DISTANCE  OF  PEAK  PRESSURE  FROM 
NOSE  OF  UNIT-CHORD  HYDROFOIL 
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FIGURE  31,  SECOND  OESIGN PROCEDURE,  PITCHING 
MOMENT  PERFORMANCE  MAP 

CAVITATION  NUMBER,  K = .l 

, T«  CAVITY  THICKNESS  AT  TRAILING 

EDGE  OF  PROFILE 
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Figure  32  - Comparison  of  Performance  Obtained  from  First  and  Second  Design 
Procedures  at  Two  Design  Lift  Coefficients,  K a 0 
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Figure  33  - Comparison  of  Performance  Obtained  from  First  and  Second  Design 
Procedures  at  Two  Design  Lift  Coefficients,  K«.05 
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Figure  34  - Comparison  of  Performance  Obtained  from  First  and  Second 
Design  Procedures,  K=>.10 
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Figure  35  - Cavity  and  Wetted  Surface  Geometry 

Elliptic  Pressure  Distribution,  Peak  Pressure  Location  s 
Second  Design  Procedure 
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Figure  36  - Cavity  and  Wetted  Surface  Geometry 

Elliptic  Pressure  Distribution,  Peak  Pressure  Location 
Second  Design  Procedure 


Figure  37  - Cavity  and  Wetted  Surface  Geometry 

Elliptic  Pressure  Distribution,  Peak  Pressure  Location 
Second  Design  Procedure 


Figure  38  - Cavity  and  Wetted  Surface  Geometry 

Elliptic  Pressure  Distribution,  Peak  Pressure  Location 
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FIGURE-40  . RLPHR  CRITICPL  VERSUS  CfiVITfiTION  NUMBER  * SECOND  DESIGN  PROCEDURE 
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FIFURE-41  LIFT  TO  ORRG  RRT10  VERSUS  CflVITRTION  NUMBER 
FIRST  DESIGN  PROCEDURE.  
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FIGURE-45  LIFT  TO  DRAG  RATIO  VERSUS  CRVITATION  NUMBER 
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FIGURE-47  ORRG  COEFFICIENT  VERSUS  CRVITRT I0N  NUMBER 
SECOND  DESIGN  PROCEDURE. 
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APPENDIX  A 


Tables  of  Performance  Data  - First  Design  Procedure 
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Symbols  in  the  following  tables  are  defined  as  follows: 


CAY 

LAMDA 

CL 


TEE 


SS 

LOD 

ECAP, 

CN 

CD 

ALPHA 

X 


= Cavitation  number,  K.  See  Equation  (1)  of  text. 

= Double  ellipse  pressure  distribution  parameter  defined 
in  Figure  2 of  text. 

= Section  lift  coefficient.  See  Equation  (3)  of  text. 

5 Cavity  thickness  parameter,  T,  illustrated  in 
Figure  1 of  text. 

= Location  of  peak  pressure,  measured  from  nose  of  profile. 

See  Figure  1 of  text. 

= Lift-to-Drag  ratio  of  profile. 

U,  DINT,  CO:  These  are  characteristic  parameters  appearing  in  the 

calculations  which  are  described  in  detail  in  Ref.  2. 

= Pitching  moment  coefficient.  See  Equation  (5)  of  text. 

= Drag  coefficient.  See  Equation  (4)  of  text. 

S Angle  of  attack  in  degrees.  Sec  Figure  1 of  text. 

£ Center  of  pressure  location,  fraction  of  chord. 
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APPENDIX  B 


Tables  of  Performance  Data  - Second  Design  Procedure 


The  symbols  in  this  table  are  the  same  as  those  of  Appendix  A 
except  for  ACAP  and  EM.  These  are  defined  in  Reference  2 as: 


ACAP  = A,  the  strength  of  the  leading  edge  singularity  due  to 
the  increase  of  the  design  attack  angle  above  the  value 
for  shockless  entry. 

EM  = m,  the  fraction  of  which  is  carried  by  the  prescribed 
pressure  distribution.  The  fraction  of  CL  carried  by  the 
angle  of  attack  term  is  1-M. 

ALPHA  = Angle  of  Attack  in  radians  (note  that  in  Appendix  A this 
quantity  is  expressed  in  degrees) . 
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APPENDIX  C 


Pressure  Distributions  for  Fully  Cavltating 
Hydrofoil  Profiles  Which  Use  the. 

Stratford  Pressure  Rise 

We  have  seen  that  fully  cavltating  hydrofoils  have  the  best  L/D 
if  the  center  of  pressure  on  the  wetted  surface  is  as  near  to  the  pro- 
file nose  as  possible  when  cavity  thickness  constraints  are  specifi- 
cally accounted  for  in  the  design.  Even  so,  there  may  be  some  situa- 
tions when  tail  loading  may  be  desired.  Then  it  is  necessary  to  avoid 
turbulent  separation  and  to  make  certain  that  it  will  not  occur  at  the 
very  beginning  of  the  design  process. 

In  particular,  if  one  specifies  the  cavitation  number,  the  Reynolds 
number,  and  the  profile  lift  coefficient,  he  needs  to  know  how  close  to 
the  trailing  edge  he  can  locate  the  center  of  pressure  and  still  main- 
tain a prescribed  value  of  the  cavity  thickness.  In  order  to  move  the 
center  of  pressure  to  the  rear,  one  needs  to  specify  adverse  pressure 
gradients  on  the  wetted  surface  which  become  steeper  and  steeper  as  the 
center  of  pressure  moves  aft.  At  the  same  time,  turbulent  separation 
must  be  avoided.  It  is  the  avoidance  of  turbulent  separation  along  with 
other  design  considerations  noted  above  which  has  led  to  the  following 
considerations.  For  nose-loaded  profiles  such  considerations  are  not  needed. 

Fortunately,  the  problem  of  turbulent-separation  avoidance  was  solved 
by  B.  S.  Stratford  in  1959  (12).  In  two  papers  in  the  Journal  of  Fluid 
Mechanics  he  has  given  a solution  for  the  form  of  pressure  rise  in  a tur- 
bulent boundary  layer  which  is  just  on  the  verge  of  separating  during  the 
whole  time  that  a typical  boundary-layer  fluid  particle  experiences  a 
ileing  pressure.  Although  such  a boundary  layer  flow  is  a very  special 
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one,  it  obviously  has  great  technical  interest,  because  it  is  not  possible 
to  prescribe  steeper  pressure  increases  without  producing  turbulent  sepa- 
ration. Of  equal  importance  for  applications,  Stratford  was  able  to  pro- 
duce such  a flow  experimentally  and  to  measure  the  Important  boundary- 
layer  flow  quantities. 

Since  Stratford’s  pioneering  efforts,  his  results  have  been  applied 
in  this  country,  principally  by  engineers  at  Douglas  [13,  14,  15,  16],  t0 
the  design  of  very  efficient  high-lift  airfoils.  Moreover,  they  have 
verified  the  correctness  of  their  methods  in  a series  of  wind  tunnel  tests 
which  showed  excellent  agreement  between  calculated  and  measured  perfor- 
mance. On  the  basis  of  these  results,  there  is  every  reason  to  believe 
that  the  work  of  Stratford  can  be  applied  to  the  hydrodynamic  design  of 
certain  supercavitating  hydrofoils. 

On  the  other  hand,  the  design  considerations  appropriate  to  cavity 
flows  are  somewhat  different  from  those  for  airfoils.  Therefore,  it  is 
necessary  for  us  to  adapt  some  of  the  methods  of  Stratford  and  of  Liebeck 
to  the  present  inverse  design  prodedure  [2].  In  this  note,  we  consider 
only  profiles  designed  in  accordance  with  the  first  profile  design  proce- 
dure. The  use  of  the  Stratford  recovery  for  the  second  design  procedure 
will  differ  from  that  given  below.  We  will  not  consider  it  hero  because 
of  its  complexity  and  its  limited  usefullncss. 

Prescribed  Pressure  Dist r i b u tlon 

Figure  C-l  shows 
schematically  the  kind 
of  pressure  distribu- 
tion under  study.  The 
wetted  surface  of  the 


Figure  C-l 
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hydrofoil  is  of  unit  chord,  and  it  lies  on  the  X-axis  as  shown.  The 
dotted  contour  illustrates  the  upper  surface  of  the  cavity  which  ori- 
ginates at  the  profile  nose.  The  lower  surface  of  the  cavity  leaves  the 
profile  trailing  edge  at  X = 1.  The  cavitation  number  h is 


^ > 

1/2P  V 

where  P^  and  are  the  steady  upstream  flow  conditions.  The  cavity 

pressure  P is  the  lowest  pressure  in  the  flow  field,  and  p is  the  liquid 
c 

2 

density.  The  disturbance  pressure  in  Ref.  2 is  P(X)  - P + p(X).J/2p  U , 

c 00 

where  P(X)  is  the  static  pressure  on  the  foil  and  the  distrubance  pressure 
is  p(X) . The  pressure  coefficient  Cp  is 

*00  - Pro 

C s — « -K  + p(X)  . 

P 1/2P  uw 

Note  that  p(X)  =0  on  the  cavity.  Moreover,  the  velocity  magnitude  on  the 
cavity  is 

q - U / 1 + X 

and,  if  all  disturbance  velocities  (u,v)  are  normalized  with  respect  to 
this  cavity  speed,  we  have 

U(X)  - qc  [1  + u(X) ] . 

and,  to  terms  of  first  order, 

o 

p(X)  - -2  (1+K.)  u(X)  . 


1 

[ 

< 

t 
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We  can  now  discuss  the  pressure  distribution  shown  above  in  terms 


of  this  notation.  Note  that  this  diagram  shows  three  regions: 

(i)  from  X = 0 to  X = Xq  at  A,  where  p(XQ)  = Xh,  X « 1. 
Two  shapes  for  this  region  are  considered; 


P(X) 


b* 

o 


or 


p(X)  = Xh‘ 


X 


linear  rise. 


parabolic  rise. 


(C-l) 


This  region  of  small  positive  perturbation  pressure  is  required  to  give 
the  profiles  resulting  from  the  design  process  some  tolerance  to  angle  of 
attack  variations  without  wetted  surface  cavitation  occurring. 

(ii)  from  A to  B [Xq  S X S (1-S)],  "Stratford  Recovery"  region. 

This  is  considered  more  fully  below. 

(iii)  from  B to  C [(1-S)  SX  < 1].  This  region  of  falling  pressure 
will  have  an  elliptic  contour; 


p(X)  - § /(l-X)  (2S-140O 


(C-2) 


This  particular  shape  has  been  chosen  to  have  zero  slope  at  B,  where 
p «*  h,  and  to  behave  like  p - /l-X  near  the  trailing  edge.  This  latter 
condition  characterizes  the  trailing  edge  condition  of  the  flows  studied 
in  [2]. 

Canonical  Pressure  Distribution: 

For  boundary  layer  calculations 
in  the  Stratford  region,  the  so-called 
canonical  pressure  distribution  is 
most  useful.  It  is  based  on  the  pressure 

Figure  C-2 
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at  p0  = P<V>  just  st  the  start  of  the  pressure  rise.  Note  that 
this  definition  varies  somewhat  from  airfoil  [13,  14,  15,  16]' usage  in 
which  Po  is  taken  to  be  the  lowest  pressure  on  the  foil.  Our  definition 
is  consistent  with  the  requirement  that,  except  for  the  leading  and 
trailing  edges,  the  wetted  surface  of  the  profile  must  have  a positive 
perturbation  pressure.  The  canonical  pressure  distribution  is  defined 
as  follows: 

P(X)  - P 

C = — - . 

P 1/2  U 2 
o 


We  can  now  translate  this  quantity  into  those  which  we  will  need  for 
hydrofoil  design.  Thus, 


C 

P 


£3 


P(X)  - P 
c 

1/2P  U*>2 


+ 


l/2p  \}J 


[pOO  - pj 


> 


where  pQ  is  the  perturbation  pressure  at  X ° XQ. 
Moreover , 


Uq2  qc2[l  + u(XQ))2  (1  + K)  [1  - po/2(l  + K))2 


■ ^ 2 “ • 

1 + K - p + p /4(1  + K) 
ro  ro 

2 

Po 

But  1 > Po  > and  we  can  neglect  the  second  order  term  compared 

to  1,  Moreover,  in  the  present  case  pQ  « Xh,  So,  the  final  result  which 
we  require  to  translate  specifications  on  perturbation  pressure  p(X)  into 
canonical  pressure  distributions,  C^(X),  is 
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Cp(X) 


p(X)  - Xh 
1 + K - Xh 


(C-3) 


J 


The  perturbation  pressure  is  needed  to  calculate  the  lift  coeffi- 
cient, C^,  and  the  canonical  pressure  coefficient  is  needed  for  boundary 
layer  calculations. 

Now  the  boundary  layer  flow  state  at  the  start  of  the  Stratford 
recovery  is  characterized  by  the  momentum  thickness  at  point  A [10,  11]. 
Moreover,  Stratford's  experiments  and  his  theory  pertain  to  a pressure 
recovery  region  which  is  preceded  by  flat  region  of  xq.  In  the  pressure 
recovery  region  normalized  coordinates  are  x/xq  with  origin  at  x = 0. 

To  start  the  pressure  rise,  Stratford  gives  formulae,  based  on  matching 
the  momentum  thickness  at  x/xq  a 1,  to  account  for  other  than  flat  pres- 
sure distributions  ahead  of  Xq  for  both  laminar  and  turbulent  boundary 
layers.  In  the  present  application,  we  will  assume  that  the  boundary 
layer  is  turbulent  from  the  nose  of  the  foil  to  the  point  Xq.  Then,  the 
"effective"  flat  region  xq  is 


rx 

0 

' ' 

3 

rx 

0 

u 

J Y a 

[2(1  + K)  - p (X) 

r\ 

u 

l 0 J 

QA  a 

* 

[2(1  + K)  - po 

A 

(C-A) 


For  the  linear  pressure  distribution,  we  get 

i2 


x “ 
o 


V f,  , 3 , Xh  5 Xh  Y 

o [_  4 1 + K 8 (l  + KJ  * ’ * J * 


and  for  the  parabolic  nose  pressure,  we  have 

>2 


x <=»  X 

O 0 


fi  + 3 Qh  V 1 

_ 2 1 + K 2 [1  + k)  J ‘ 
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In  either  of  these  cases,  we  can  write 


x ® X — 

O O T 


(C-5) 


The  quantity  T is  defined  by  one  or  the  other  of  the  above  functions 
of  Xh/(1  + K). 

The  Stratford  Recovery 

The  chief  result  of  Stratford’s  investigation  is  contained  in  his 
canonical  pressure  distribution  [12],  it  is 


t 


0.49  {Re  1/5  [(  — )1/5  - 1]}1/3,  C < 4 
o x p 7 

o 

(C-6) 

1 - , C > y 

\/  (x/xQ)+b  ^ 


where  Re  « 
o 


U<”Xo  /l'+  K 
V T 


and  where  the  constants  a and  b are  chosen  to  match  C and  dC  /dx  when 
_ P P 

0^  » 4/7.  Now,  the  above  formulae  are  written  for  a hydrofoil  of  unit 
chord.  Naturally,  we  will  be  dealing  with  foils  of  chord  c.  All  we  need 
to  relate  Re^  for  this  case  to  the  Reynolds  number  based  on  chord  c is 
to  put 


fU  Cl 

CO 

Re  **  — — 
o v 


U O',  X A + K 


*»  Re 


x AT+k 

o 


(C-7) 


where  now,  Re  and  Re  are  based  on  actual  distances. 

o 

Coordinates  and  Lift  Coefficient 


As  we  have  seen,  there  are  two  sets  of  coordinates:  those  which 
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measure  actual  distances  along  the  unit  chord  of  the  profile,  designated 
by  X,  Xq  etc.,  and  those  which  relate  to  the  canonical  pressure  distribu- 
tion, designated  by  x,  xq  etc.  We  must  determine  the  relationships  be- 
tween these  coordinates  in  order  to  continue  the  analysis.  We  have  al- 
ready obtained  one  result  in  the  form  of  the  Equation  (C-5)  which  is 


X - x T . 
o o 

However,  this  relationship  holds  only  at  the  point  A on  the  foil.  That 
is,  when  X = Xq  and  x = xq.  But  we  know  that  in  going  from  A to  B,  we 
follow  the  Stratford  Recovery  and,  in  particular,  when  we  are  at  B, 

X = (1  - S) . Moreover,  at  B we  know  that  p(X)  = h,  and  so  the  canonical 
pressure  coefficient  from  Equation  (3)  is 


C 

P 


i).h- 
1 + K - Xh 


(C-8) 


This  value  of  can  be  put  into  the  Stratford  Recovery  and  the  cor- 
responding value  of  x/xq  found.  Denote  this  value  of  x/xq,  corresponding 
to  X ° 1 - S,  by  Z.  Then,  we  find  that 


Z » 


1 + 


<keo) 


r (i  - x>h  t 

I7T  |_.49(1  + K - XhTJ 


(1  - X)h  4 

1 + K - Mi  7 


r (i  - xyr 

x K Xl- 


• i 


+ K - Xh 


> 


4 

7 


(C-9) 


We  now  have  two  pairs  of  corresponding  points; 

x . X . . . „ x „ X 1-S 

at  A — * 1,  rr*  •»  1;  and  at  B * — * Z,  — » . 

x * X xt  ‘X  x 


Evidentially,  as  one  traverses  the  distance  from  A to  B on  the  wetted 
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surface  x goes  from  x to  Zx  , Therefore,  the  distance  is  (Z-l)x  = 

oo  ’ o 

(Z-1)Xo/t.  But,  the  profile  has  unit  chord.  Hence,  the  sum  of  the 
distances  from  0 to  A,  A to  B and  B to  C is 

1 - X + (Z-l)X  /t  + S 
o o 

or 

(l-S)x  - <T  + Z - l)Xo  (C-10) 

This  is  a fundamental  geometrical  constraint  on  the  variables  of  the 
problem  which  we  will  need  later. 

In  order  to  consider  further  the  relationship  between  the  boundary 
layer  coordinate  x and  the  hydrofoil  coordinate  X,  let  us  take  a point 

Q in  the  Stratford  Recovery  region  between  A and  B.  Suppose  Q has  the 

foil  coordinate  X and  the  boundary  layer  coordinate  x.  Then,  it  must 

follow  that 

X - X a x - x 
o o 

Now  we  can  divide  this  relationship  by 


Next,  we  can  apply  Eq.  (C~5)  to  this  result  to  obtain 

“x~  " + 1 “ T * (C-ll) 

O O 

Equation  (C-5)  relates  the  boundary  layer  variable  z *■  x/xq,  appearing 
in  the  Stratford  Recovery  equations,  to  the  geometric  coordinate  X on 
the  profile. 
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It  is  the  coordinate  X,  along  with  the  perturbation  pressure  p(X), 
which  allows  us  to  calculate  the  lift  coefficient.  Moreover,  since  the 
quantity  is  determined  by  Eqs.  (C-l)  or  (C-2) , we  can  use  it  to  relate 
the  boundary  layer  coordinate  in  the  Stratford  region  to  the  geometric 
coordinate  X without  ambiguity.  In  particular,  we  can  write 


C 


L 


fl 

p(X)dX 

-0 


Then,  with  the  help  of  Eqs.  (C-l),  (C-2),  and  (C-3) , this  formula  becomes 


CL  ° Xh[(l-S)  ~ | X0)  + ^ * (1  + K - Ah)  — 


T4-Z-1 

C (—  + 1 - T)  d<^-) 
p xo  xo 


(C-12) 

for  the  linear  nose  pressure  distribution  and  for  the  parabolic  case, 


irsh 


X 


\ - AMI  - S - f Xo]  + -fi  + (1+K-Ah)  - 


T+Z-l 

Vr + 1 - a<r> 


o 
fr-i 


Note  that  Eq,  (C-ll)  and  the  considerations  leading  to  Eq.  (C-10)  hove 
been  used  to  put  the  integrals  in  Eqs.  (C-12)  and  (C-13)  in  a form  suit- 
able for  direct  use  of  the  Stratford  Recovery  as  specified  by  Eqs.  (C-6) 
Further  Details  of  the  Solution 

The  preceding  sections  provide  all  essential  ingredients  necessary 

for  the  determination  of  the  prescribed  pressure  distribution  on  the 
profile.  However,  in  order  to  implement  these  results,  a number  of  sub- 
sidiary calculations  oust  bo  completed.  For  example,  we  have  ‘noted  that 
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the  Stratford  Recovery  requires  that  a joining  of  profile  shapes  occurs 

at  C — 4/7 . At  this  pressure,  C (z)  and  dC  / dz  are  to  be  continuous. 

P P P 

These  two  conditions  permit  us  to  calculate  the  constants  a and  b of  Eqs. 
(C-6) . At  the  point  of  joining,  let  us  put 


z = z for  C = -7  . 

o P 7 


Then,  from  the  first  of  Eqs.  (C-6),  we  have, 


z * 
0 


1.586 


1/5 


Then,  it  follows  that 


1l?M  , */5  , 1/5  _ .2/3 

175  o a ' 


b =■ 


Re 


, ,,,  A/5  , 1/5  , .2/3 

6.563  z (z  - 1) 
o o 


Re 


1/5 


Now,  we  can  continue  with  the  determination  of  the  matching  of  the 
Stratford  Recovery  to  other  prescribed  conditions  on  the  pressure  distri- 
bution. We  require  that  the  lift  coefficient  be  specified  in  advance. 
However,  the  point  Xq  is  as  yet  undetermined,  and  the  value  of  h,  the 
maximum  pressure  on  the  wetted  surface,  is  also  unknown.  We  know  that  the 
Stratford  Recovery  allows  for  the  steepest  possible  pressure  rise  without 
turbulent  separation.  Thus,  we  can  specify  a family  of  pressure  distri- 
butions, depending  on  the  parameter  S.  Members  of  this  family  will  all 
have  the  same  values  of  C^»  K and  Reynolds  number. 

The  values  X and  h must  be  determined  by  trial  and  error.  Wc  can 
o 

use  the  basic  geometrical  constraint  and  either  Eqs.  (C-12)  or  (013) 
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for  to  guide  the  iteration  in  a Newton-Raphson  procedure.  Thus,  we  can 
write  Eq.  (C-10)  as 

G <XQ,h)  = “--l-Z  + l»0  , 

° o 

and  Eq.  (C-.12) , for  example,  as 

C.  (X  ,h)  » X[1  - S - ~ X ] + ttSh/4  + I 

L 0 l O 


Then,  we  suppose  that,  in  the  process  of  iteration,  we  find  at  some 

step  a difference  in  the  calculated  value  of  CT  from  that  specified. 

Li 


Let  this  difference  be  AC, 


Then  we  have 


and 


8C  3C 

aiT  Ah  + aif" 


AC, 


f Ah+|f  AX  « 0 
3h  dX  o 
£> 


Evidentally , one  can  solve  for  Ah  and  AX  and  use  them  to  determine  the 

o 

next  value  of  h and  Xo  in  the  iteration.  This  is  a standard  procedure. 

If.  needs  no  further  discussion  here.  The  calculation  of  all  the  partial 
derivatives  from  the  formulae  of  the  preceding  sections  is  rather  tedious. 
We  shall,  therefore,  proceed  to  the  discussion  of  results. 

These  results  are  contained  in  Figure  (C-3)  which  shows  the 
variation  of  L/l>  for  a range  of  peak  pressure  locations  at  zero  cavita- 
tion number.  For  each  value  of  cavity  thickness  T three  values  of  are 
used.  Note  that  the  trends  show  that  the  nose  loading,  corresponding  to 
large  values  of  S,  lead  to  more  favorable  L/D  values  than  the  tail -loaded 
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profile  obtained  at  small  values  of  S.  As  we  have  already  explained  the 
chief  value  of  the  Stratford  recovery  is  for  tail-loaded  profiles  which 
are  not  as  efficient  as  those  having,  nose  loading.  Therefore  we  have 
discarded  this  kind  of  pressure  distribution  in  favor  of  the  simple 
double  elliptical  form  which  has  been  used  in  the  body  of  this  study. 
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APPENDIX  D 


..A 


Fig*  1 Schamolk  diagram  of  cavity  llow  pool  a hydrofoil 


tho  r-nxis.  Tho  chord  of  Iho  profile  is  inclined  At  An  angle  of 
attack  « with  rcs|W>ct  to  tho  undisturbed  flow.  Tito  shape  of  tho 
welled  Btirfnco,  measured  from  tho  r-Avis,  is  defined  by  the  func- 
tion/(x).  Tho  cavity  extends  back  over  tho  upper  .surface  of  tho 
foil.  Generally,  tho  statin  pressure  of  tho  gas  or  vapor  iiivdo  tho 
cavity  I\  will  bo  less  than  the  free  stream  static  pressure  and 
tho  cavity  will  close  at  a finite  dislanro  of  abseiMsa  f,  hehitid  tho 
hydrofoil  noao.  Whon  I > 1 Iho  How  is  called  fully  cavitating. 
Partial  cavitation  is  said  to  exist  when  tho  cavity  rluym  on  tho 
uppor  mtrfaeo  of  tho  profdo.  Otdy  fully  caviUiting  flmvs  uro  cua- 
sidored  in  this  paper.  In  such  flows  tho  How  state  is  characterised 
by  Uto  CAvithliou  number, 


<1) 


Hydrofoils  at  Zero  Cavitation  Number  and  Related  Airfoils* 

Consider  the  steady  rectilinear  flow  of  a liquid  of  density  p past 
4 fully  cn vita  ling  hydrofoil  of  unit  chord  in  two  dimensions.  As 
'Uustrated  in  Fig.  1,  tho  profdo  noeo  is  at  tho  origin  of  coordinates 
\d  tho  freo  stream  velocity  luu  magnitude  U and  flows  parallel  to 

• Numlirf*  in  hrnrkcu  designslo  lloferoiircs  si  end  of  pnpor. 

•Tim  resell*  o(  ibis  Mvlttm  are  duo  to  T.  V.  Wu,  who  hns  kindly 
pertained  them  to  Ik:  used  Imro.  Professor  Wu’s  helpful  diseuwton  of 
many  Aspect*  of  tht*  p»i>er  is  yr *lofully  seknowlodgod. 

Contributed  tiy  tbs  llydrstdio  DivUInn  and  prevmltd  at  ths  Winter 
Annual  Meeting,  PluUdcIphiA.  Pa.,  N’ovombor  IT— 33,  1003,  of  Ttts 
Amc.nican  Hocic.rr  or  Mkcuanioal  limn  ns***.  Manuscript  rs- 
otivsd  at  Ati.Mtt  llssdqusnsr*.  May  13, 1064, 
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In  this  section  we  shall  consider  thnjo  rpcri.nl  flows  for  which 
K - 0.  In  this  limiting  caxo  tho  envity  length  f lKtroir.es  mlunlo 
snd  the  flow  sjxtod  oo  tho  cavity  surfAco  is  ctpial  to  the  free  stream 
speed,  U, 

Let  the  flow  st  sny  point  in  liic  x,  y-platte  bo  given  by  tho 
velocity  vector  q(r,  p),  Next,  as  is  customary,  wo  introduce 
tho  2 and  p components  of  the  disturbance  velocity  mused  by  tho 
foil.  These  are  u (r,  p)  and  i<x,  y),  rcs|veclivoly.  Tho  perlorU-i- 
tion  velocities*  sra  reiated  to  tho  vector  q by 


q «•  {/|1  + «i,  v|. 


(2) 


In  accordance  with  tho  liitcarijuilhm  scheme  wo  must  have  u,  r « 
1.  Wo  shall  also  Sasunia  that  tho  cavity  surface*  and  the  actual 
surface  of  tho  hydrofoil  lio  chwo  enough  to  the  X-axis  «,i  that  in  iho 
x-plano  (s  ••  C 4-  tp)  Uio  boundary  conditions  on  tho  complex 
velocity, 


ts(s)  » u - it. 


<3> 


occcmxi  i s s < / 641 
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Flg,  3 Doundory  condition,  In  lh»  phyilcal  plan*  U-plans)  ond  con*, 
•ponding  boundaries  in  1h»  Iramiwmod  plan*  (vplons)  It  !«•  c*vll*» 
ilonnvmbw 


cnn  lot applied  on  llio  x-nxis,  Tho  linonriml  boundary  conditions 
/or  tho  analytic  function  u{t)  nro  shown  in  Fig.  2.  Theso  condi* 
lions  can  bosuinmnriicd  ns  follows: 


u(x,  0+)  - 0 for  x > 0,  . ' 

v(x,  0- ) - /'(*)  for  0 <,  x £ 1, 
u(x,0-)  -0  for  x > 1, 
u - o - 0 as  |*|  -*  m. 

liy  means  of  tho  conformal  transformation, 

f - £ + it  - 


(4) 


(6) 


tho  wholo  r-plano  is  mapped  into  llio  lower  half  o*  tho  s-plano 
with  tho  solid  surface  lyinj?  on  tho  real  y-axis  between  £ - 0 ond 
£ <•  l.  Tho  cavity  boundaries  lio  along  tho  remainder  of  tho  real 
axis.  The  r-plmto  and  tho  transformed  boundary  conditions  are 
also  illustrated  in  Fig.  2.  It  is  to  bo  noted  that  wo  lutvo  chosen  to 
to  be  invariant  under  tho  mapping,  that  is 


to*(P)  «•  — to*(v), 


V m f(|)* 

Then  from  s^ualion  (0)  it  follows  that 


lienee,  tho  shape  of  the  related  airfoil  is  given  by 
61?  t stceuii 
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Fig.  3 0»*m«lrl<  pwgmtlMI  far  a clrculor-*rc  hydrofoil 


ett) 


- 

Jo  2‘ 


‘ 1 W'>,u 

<U 


• 1 (0) 


«<*)  - u(x,  p)  - iv(Xf  p)  - wfy')  - w*(y) 

- u*(£, »?)  - «>*(£.  «f)<  <Q> 

Wo  ran  obtain  tho  relaletl  airfoil  formulas  by  noting  that  on 
part  of  tho  real  p-axis  u — u*  “ 0.  Theroforo  w*(y)  can  bo  eon* 
turned  analytically  into  tho  upper  half  of  tho  y-plano  by  tho 
Itiemann-Schwnrls  principle  of  reflection  (0).  For  tho  proeeat 
flow  wo  can  write 


for0<  £ < 1. 

For  oxnntplo,  if  we  tnko 


/(*)  ••  o*  + bx*  ex*  + . 
we  find  that 


0 £ 1, 


»<£>  ~ «£  + £•  + * 0 S £ < 1. 


Kviilenlly,  an  inclined  flat  plate,  detini’il  by  f(t)  " «/,  is  mapped 
into  itself.  In  tho  general  raso  the  shape  of  tho  rnluled  airfoil  lind 
ovon  its  angle  of  attack  will  bo  different  from  thoou  of  lltc  original 
hydrofoil. 

Corresponding  to  tho  disturbance  velocity,  wo  shall  define  tho 
perturbation  prosauro  coofikiont  pfx,  p)  by 


. P(r,  v)  - 1\ 

p(x,  p)  ( , 


00) 


- pf/* 
2 


where  /'(x,  p)  denntm  tho  static  prenaiirn  lit  any  pmnt  in  llio  llntv. 
Uso  of  tho  linearised  nnrooulli  cipiatiim  allows  that  jifx,  p)  •• 
— 2u(x,  y).  Tho  lift,  drag,  and  moment*  caolhcknls  nro 


/ CL\  ri  ( 1 \ 

( Co  ) <-  I (/'(*) ) pfx,  0-)<ix 

\cj  Jo  \-x/ 

•»  — 2 J'  ^v'(x,  0-)j  u(r,  0- >ix 

* /o'  (*•<*•  «,->)  ' ( i i? 


llio  expression  for  Ca  can  bo  simplified  by  means  of  tho  I'oituttm 
integral  formula  (7J  • ’ . 


(7) 


so  that  u*  is  odd  and  u*  is  even  in  » ).  After  tho  analytic  con* 
Imuation  the  flow  in  the  e-plnno  will  be  recognised  as  tho  flow 
)Kist  an  airfoil  of  tern  thickness  with  ui*  being  tho  diaturbaneo 
created  by  tho  airfoil  in  an  otherwise  uniform  stream  of  velocity 
7*  «•  {/.  As  a emwctpicnco  of  tho  invar ianro  indicated  by  orpin* 
tion  (0),  the  slm|m  of  tho  transformed  solid  boundary  in  tho 
p.plnno,  determined  by  v*(f,  0),  will  generally  bo  different  from 
Uto  welled  surface  sltapo  in  thea-plane.  Lot  tbs  shop*  of  Ute  now 
related  airfoil  bo  denoted  by 


£ 


dl, 


where  P.V.  alaadt  for  Uto  Cauchy  prineiplo  valuo  of  Uie  integral. 
Thus, 


' !(£-<)  t\ii> 


-TM-  />“■' 

- “ |J'  u»(£,  0->d<}*  + <JJ  v»«,0-)-V,o-)fdi 


(8)  Hen  co 


~~  - *x,  0-)  - i*(£,  0*)  - for  0 <*-£•<  1. 


(W> 


• Momenu  tro  rountol  |k»Mtivo  if  lluiy  Irutl  Iti  mrtrrtM*  tin*  .*»t»;lo 

of  MUrk*  Koto  tkoi  iUo  tuotuout  is  Uko«  rIh»ui  Uto  o*  Uh> 

pro&lo. 
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n»  tlio  ill  her  Imtul,  the  lift,  moment,  ami  tho  third  .muncul  of  tho 
■\elnlcd  airfoil  ore 

. ^ Q (;**((.  0-)  - vU  o+y,  t 

u*((,  0— (13) 

Therefore,  by  comparison  with  equations  (11)  and  (12)  we  find 

C,.  - (V,  <V-  P,  - <V*/2ir.  (14) 

The  loon  (ion  of  the  renter  of  pressure,  oxpnxwod  iui  ft  fraction  of 
the  chord  from  (ho  hydrofoil  nose,  is 

C.l\  - -Cw/C’i,  - -C,„*/C„*.  (Ho) 

These  equivalence  relationships  between  tho  hyd’  ofoil  and  tho  re- 
lated airfoil  are  duo  to  Tulin  and  Burkart  (1,  2). 

From  tho  direct  problem  of  thin  airfoil  theory  18),  it  is  known 
that  CL* i Cm*,  and  C,h*  can  bo  ex>  ■ ceeod  in  terms  of  tho  Fourior 
coefficients  of  v*(£,0=fc).  Write 

m 

-v*(f,  0±)  - ~ 2 °» 000  nff*  <16) 

4 n- 1 


hyihiifiiil.  Nest,  In  show  ( In*  elleilit  nf  l aiiilu  r,  we  will  npp 
(.henry  lo  circular  air  liydiofodu.  Fittollv,  v.e  mil  luoe.iiii 
(i|ilimiMO  design  pruiilein  fur  fully  rnvit.'iling  hvdtnfnds  at  A 
Exempts  t.  ri«l-PloU  Hyd.sfolt.  Fora  llal  plaUi  hydrofoil 
angle  of  attack,  a,  wo  Itnvo/fr)  •»  - it/  and 

«*(£,<>-)  •*  Kr,0-)  « 

From  equations  ( 1 1 ) it  follows  that 

04  «•  2 or,  <«,  - 0,  » « 1,  2, ... . 

Ilone.o,  from  equations  (M)  and  (10)  wo  ohla'm 

Cl  - J «,  C„  - ' «*,  . 

Cm  - - '!*  « mid  0.1*.  " jjl' 

These  results  hhow  the  total  force  on  the  fully  eavil-'iting  flal 
to  he  normal  to  tho  plain,  as  I1-  should  lie.  Also  O.Ih,  the  :d 
of  tho  center  of  pressure,  is  */»  rhoril  from  I lie  leading 
Thcso  values  for  the  coelficienls  will  be  rerngnixed  ns  a line 
form  of  KirchofT’s  exact  solution  [91  for  small  values  of  2. 

Example  I.  Clrtular-Are  Hydrofoil.  The  circulur-IiT.  ] 
which  subtends  a total  anglo  of  2y  at  tlin  center  of  the  > ire 
which  ia  at  an  anglo  of  attack  a,  is  defined  for  small  values  c 


where 


f - (1  + cos  d)A  ~T<$<  V,  ..  • 


and  where 


a.  „ _!  P '„•({,  Odb)  coe  nff  69,  (16a) 

* JO  s ' 

I 

It  is  also  ehown  that  . . 


0 


p*(£, 0-)  w-  -2w*(£,0-)  - -astnn- 


— 2 ^ a.  elnnd,  — r<0<O,  (106) 

* • l 


and  hence,  from  equation  (13)  that 


Ci*  - *(««  + Oi),  c7 u*  " ~ (a*  + 2*  + o«)j 


ai  1 


Ciu*  •*  t*“  <4  F 7(<j|  + a<)  + 3aj  + — 


(10) 


/(i)  - -02+72(1  - x). 


Then  wo  have 


*(f,0-)  - t<x,0-)  - 7 - « - 272  - 7 “ <*  - ,2y{ 


-(?-)- 


7 c«s  0 — 


| 4 


wlmro  9 ie  definod  In  equnlinn  (15).  From  equation  (11) 
tain  . . . 


a,  2o  - oi  - 7i  o»  - 7/4, 


o,  - n»  0 

From  equations  (14)  and  (10)  it  follow*  tlmt 


and 


Thcrcfiire,  from  cquntWnii  (14)  ami  (10),  tho  hydrodynamic  mao. 
tie  is  on  tho  hydrofoil  sro  determined  by  quadrature*. 

It  is  important  In  note  that  for  vaporous  cavitation  the  rnvity 
pressure  rt  -ist  lie  the  lowest  pressure  in  tho  flow.  Thcrcfnre,  <m 
tho  welled  surface  of  tho  foil,  the  disturbance  pressure,  p(x,  0-) 
•»  ;>•((,  0- ),  most  be  nnnnrgativo.  Accordingly,  Ute  rnenieient- 
o.  must  be  such  that  tho  righl-lmnd  side  nf  equation  (156)  Is  mm- 
ncgalivc,  Otherwise  thn  llow  emdigumlinn  will  not  correspond 
to  (lint  which  we  have  pustulated.  (This  requirement  is  two 
which  should  lm  added  lo  ihnso summarised  in  equation  (4).)  In 
particular,  it  Is  orccemiry  (hat  fl»  J 0.  This  rocdirienl  dnlrr* 
tiiiues  the  strengih  nf  the  linsn  singularity . When  it  vanishes,  tho 
flow  is  said  lo  exhibit  smooth  nr  "slmekhW  entry.  If  m were  lo 
become  negative  the  cavity  would  change  sides  at  least  near  tho 
profile  leading  edge. 

Tho  preceding  development  nf  the  linearitod  theory  for  A •*  0 
pcrinils  us  to  consider  three  specific  examples.  First,  In  order  lo 
givo  a I lasts  for  subsequent  cunqiarison  w«  will  treat  tin)  flftt  pi*  to 


a,,*"l«(W+8V)/("  '■  ^)' 

Moreover  tho  requirement  Uinta*  ^ (I  implies  | lull 


if  tlio  sold l mu  it  to  lm  valid.  When  this  cnmlilion  is  sali 
see  that  at  a given  value  nf  rr,  nn  ItirroHning  pusi’.ivi* 
cnndier  always  iucrenses  (lm  lift,  drug,  end  mumcul  mu 
moves  (Ihi  cxmler  nf  pressure  further  In  Ibo  fc«r  from  (be  ( 
hwftlimu  In  parlirulnr, 


«».R,  5&.»f.+.A  r-..9 

4y  10  Oy  i \ t / oy  o 


A point  nf  Interest  Is  tho  faveralde  im-rouso  of  C ,,  duo  to  it 
7 for  fully  raviiMiug  hydrofoiU.  For  lival  «,  tho  lift'd 
Is  ft  maximum  at  1 
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M ! 


{ ) 


y!.»»  “ a, 


(25) 


giving 


40  _1 
40  a' 


(20) 


This  vnluo  of  L/D  is  a 22.5  percent  increase  over  tho  flat-plnto 
value  of  l/«. 

ExompU  3.  Optimum  Profll*  D*iijn.  From  tho  results  of  tlm  pro- 
celling  example  wo  arc  led  to  consider  tho  design  of  tho  optimum 
profdo.  To  do  this  ono  attempts  to  dotormino  that  profilo  shapo 
which,  for  prescribed  drag  coefficient  and  attack  angle  a,  will 
produco  tho  highest  lift.  From  tho  gcnoral  formulas  wo  have 


Ct  ■ — (a«  *4-  2or  + o») 


(27) 


and 


Cd  m ^ (<**  + «»)’. 


(28) 


.When  tho  value  of  a is  prescribed  wo  must  havo/(0)  — 0,  /( 1 ) ■* 
—a  which  constrains  tho  form  of  tho  Fourier  coefficients  as 
follows: 


/(l)  - -a  - v(s,  0-)dx  - 2 /;  u‘(£,  0-}£d£, 


O#  V"'  C»«  XT'' 

° “ *2  ~ 4*  4n‘  - 1 “ Y (2n  4* 


o«,.t 


l)(2n  - 3)’ 


a*  + 2 0. 
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K . . 

2 V 5*  3rr, 


and 


<33 


If  tlto  flew  configuration  is  to  remain  im  p<Mitlai.ed,  tin*  first  in 
equality  shows  that  6'/,  and  a cannot  tmtli  lie  prem-ribed  nrlii 
trnrily.  Wlion  this  first  condition  is  satisfied  the  second  inequal 
ity  gives  tho  largest  attainable  valun  of  CL.  To  proceed  with  th 
analysis  lot  us  write 


# 


tot 


and 


(34 


or  upon  using  oqualion  (16)  wo  find  after  term  by  term  integration 
that 


whero  t % 1 and  0 < i ^ 1.  For  a fixed  value  of.  6,  it  is  fmiii 
that  C/./C’/,  in  a maximum  when  t 4/3.  The  limit  S ->  1 lead 
to  fle  •=•  0,  which  corresponds  to  tho  smooth  entry  ('iitifig'-ritioi 
In  practice  it  may  he  desirable  to  retain  the  lending  edge  „ingi 
larily  sinco  smooth  entry  is  difficult  to  maintain  or  even  t 
renlir.o.  To  retain  tho  nose  singularity  we  simply  choose  6 slight! 
leas  than  unity. 

For  instanco  >f  wo  tako  < ■=*  4/3,6  *-  14/15,  then  for  engincorir 
purposes  wo  have  tho  following  optimum  values: 


■C*-f  “* 


(29)  and 


:(3i 


Moreover,  sinco  p(x,  0— ) must  bo  nonnegative  at  every  point  on 
tho  foil  wo  have  from  equation  (166),  that  <4^0,  and  that 


7ir 

Cu  - To-  «• 


~ tan  ^ -f  ^2  °*  ®‘nn  4 ? 0 for  0 <i  9 < *r.  (30) 


The  ratio  Cu/Cn,  • 7/(oa)  is  greater  than  tho  flat  plate  easo  I 
40  percent  and  it  surpasses  tho  maximum  L/D  of  the  circular  u 
by  14.4  percent.  Now  from  equations  (35)  and  (32)  wo  oUlftii 


It  is  clear  on  physical  grounds  that  if  tho  whole  underside  of  tho 
profile  Is  to  bo  wetted,  from  the  nose  to  the  trailing  edgo,  wo  must 
havot<l}  <,  0,  or  by  equation  (15)  that 


a, 


8 . 4 

oi  - •••  a and  a,  — — a. 
5 5 


(3 


(31) 


Finally,  tho  pressure  condition,  oqualion  (30),  boemnes 
«4 


_ tan  — + a,  sin  0 + at  sin  20 
2 2 


Since  Ci,  ami  Cn  depend  upon  tho  first  three  coefficients  only,  tho 
problem  of  finding  a imiquo  set  of  values  of  tho  a„*s  correopomling 
to  the  liighcst  Ct,  for  pn-scrilieil  Cn  and  «,  and  which  also  satisfies 
conditions  (30)  and  (31 ) m not  a mini. inn*  prnblmn. 

For  practical  engineering  purposes  wo  shall  confine  ourselves 
to  tho  special  class  of  profiles  for  which  <>i  •«  o«  — 0.  Then 

equations  (27)  through  (20)  can  bo  solved  for  <4,  a(l  and  <4.  Tho 
solution  is 


1 


•»  — cna  0(1  — cits  0X1  + 3 sin'  0)  > 


for  0 < 0 < r,  so  that  this  condition  In  fullillml,  Tim  enrrviqKiit 
ing  shapo  of  tho  wotted  surface  is  found  U>  bn 


/<*) 


- -ya  + a,x(l-|v'^)-“*(3-10v'x‘  + i: 


80. 

r 


BC 


0,  - — — 6 

r ! » 


r 

'nc0 


— ta, 

+ 0 a, 


and 


bCb 


<4  “ — + 4 “ 6a. 


a 


- T*(ll  -32  V*  + 10*). 

O 


0 < * < 1.  (v 


(33) 


Tlie  camber  function  of  tho  wotted  surface,  measured  from  t 
chord  is 


10a 


M*)  - /(*)  + a*  - *(1  - V*  )*. 

0 


()<*<i.  (; 


The  condition  Dial  <4  0 and  condition*  (30)  and  (31)  laad  to 
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The  last  point  to  cheek  is  to  see  whether  or  nut  there  is  hit 
feronco  between  the  upixir  surface  of  tlm  foil  and  tin*  vombor  li 
From  tho  general  solution, 


Transactions  nf  thn.  AS* 


V 


0.0 


0.4 

0.2  * ' 

y ■ 

0.1  - 


0 • 


0* 


(<0, 


(Cw), 


• •|Ja  • 0.(910 
. •IS.04 

"US0  **o.oje»4 


ftp,  4 Profit*  «nd  cavity  centaur*  at  • law-drag  hydrafall  itK>0, 
Not*  magnified  vertical  (Cat*. 

w*(xj  « u*  - W*  - ~ - ^V~7^) 

‘ mi 

+ <'£«»  (2k-  1 -2%/kk-  I)]*, 

• i 

in  which  tho  branches  arc  chosen  to  make  — l)/»  -*■  1 and 
* VKt'  — 1)  — * v m [v|  -*  »».  Wo  find  tho  uppor  and  lower’ 
boundaries  of  tho  cavity  from 

*♦(*)-  J*  v(x,0+)<i£  for  *>0, 

and 

• l . 

"*  A_(x)  - -a  + J*  v(x,0-)dg  for  1. 

When  the  values  of  a.  for  this  example  aro  put  in  tho  general 
solution  and  when  the  values  for  v(x,  0=b)  are  subs  It  luted  in  tho 
integrals  ah'  'o  tho  final  result  io 

A«(*}  - "■ ((7  - 2 y/x  - 48*  - 32**/')  y/x  + y/i 

- 7 lug  (VtT  V*  + *‘A)  + 22*  + 04*'/.  * 32**},  (39) 

M*)  - « + •—  {(7  + 2 - 48*  + 32*,A)  W-  V* 

• + 7 log  (VV*  - 1 + *'/‘)  - 22(x  - 1) 

, + 04(*'/»  - 1)  -32{*»  - 1)}.  (40) 

The  prvGlo  and  tho  up|ier  surfaco  of  tho  cavity  calculated  from 
equations  (38)  ond  (39)  aro  shown  to  scale  in  Fig,  4.  It  is  Boon 
that  tho  cavity  does  clear  tho  lower  surface  of  tho  hydrofoil. 
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